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Chapter 1 

Electromagnetic field 
fundamentals 

1.1 Introduction 

This chapter starts with a brief review of Maxwell's equations, which are the 
fundamental laws that, together with the theory of electromagnetic behavior 
of matter, explain on a macroscopic scale the properties of the electromagnetic 
field, the relationships of this field with its sources, and its interaction with 
matter. The reader is assumed to be familiar with these equations at least at an 
undergraduate level. Next, after reviewing other fundamental topics such as con- 
stitutive parameters and boundary conditions, we apply the energy-conservation 
law to a bounded volume, limited by a surface S, inside of which there exists a 
time-variable electromagnetic field. We shall see that when the energy balance 
is formulated, there appears a term representing a flow of energy carried by the 
electromagnetic field through the surface S that limits V. This term leads us 
to the definition of Poynting's vector. Similarly, when the law of conservation 
of momentum is applied to the same region, we find that the electromagnetic 
field also carries a momentum density, which can also be expressed in terms of 
Poynting's vector. 

1.2 Review of Maxwell's equations 

The general theory of electromagnetic phenomena is based on Maxwell's equa- 
tions, which constitute a set of four coupled first-order vector partial-differential 
equations relating the space and time changes of electric and magnetic fields to 
their scalar source densities (divergence) and vector source densities (curl) 1 . 

According to the Helmholtz theorem a vector field K is uniquely determined by its di- 
vergence and curl if they are given throughout the entire space and if they approach zero at 
infinity at least as 1/r™ with n > 1. A proof of this theorem is given in Appendix ?? 
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Maxwell's equations are usually formulated in differential form (i.e., as relation- 
ships between quantities at the same point in space and at the same instant in 
time) or in integral form where, at a given instant, the relations of the fields 
with their sources are considered over an extensive region of space. The two 
formulations are related by the divergence (??) and Stokes' (??) theorems. 
For stationary media 2 , Maxwell's equations in differential and integral forms 

are: 

Differential form of Maxwell's equations 

V • D{r,t) = p(f, t) (Gauss' law) (1.1a) 

V • B(f,t) = (Gauss' law for magnetic fields) (1-lb) 

VxE(r,t) = - dB g t '^ (Faraday's law) (1.1c) 

_ _ dD(ft) 
V x H(f,t) = J(f,t)-\ -^-^ (Generalized Ampere's law) (l.ld) 

Integral form of Maxwell's equations 

<j> D(r,t) ■ ds = Q T (t) (Gauss' law) (1.2a) 
<j> B(r,t) ■ ds = (Gauss' law for magnetic fields) (l-2b) 

<j> E(r, t)-dl = -J ^ ■ ds (Faraday's law) (1.2c) 

(f H(r, t) ■ dl = f (J(r, t) + ——J-^—) ■ ds (Generalized Ampere's law) 
Jr Js 

(1.2d) 

Maxwell's equations, involve only macroscopic electromagnetic fields and, 
explicitly, only macroscopic densities of free-charge, p(f,t), which are free to 
move within the medium, giving rise to the free-current densities, J(r,t). The 
effect of the macroscopic charges and current densities bound to the medium's 
molecules is implicitly included in the auxiliary magnitudes D and H which are 
related to the electric and magnetic fields, E and B by the so-called constitutive 
equations that describe the behavior of the medium (see Subsection 1.2.2). In 
general, the quantities in these equations are arbitrary functions of the position 
(r) and time 3 (t). The definitions and units of these quantities are 

E = electric field intensity (volts/meter; V mT 1 ) 



2 In a stationary medium all quantities are evaluated in a reference frame in which the ob- 
server and all the surfaces and volumes are assumed to be at rest. Maxwell's equations for 
moving media can be considered in terms of the special theory of relativity, as shown in 
chapter ??. 

3 Throughout the book, in most cases, in order to make the notation more concise, we will 
not explicitly indicate the arguments, (r, t), of the magnitudes unless we consider it convenient 
to emphasize the dependence on any of the variables. 
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B = magnetic flux density (teslas 4 or webers/square meter; T or Wb mr 2 ) 
D = electric flux density (coulombs/square meter; C m~ 2 ) 
H = magnetic field intensity (amperes/meter; A to -1 ) 
p = free electric charge density (coulombs/ cubic meter; C mr A ) 
Qt = net free charge, in coulombs (C), inside any closed surface S 
J= free electric current density (amperes/square meter A mr 2 ). 
Three of Maxwell's equations (1.1a), (1.1c), (l.ld), or their alternative inte- 
gral formulations (1.2a), (1.2c), (1.2d), are normally known by the names of the 
scientists who deduced them. For its similarity with (1.1a), equation (1.1b) is 
usually termed the Gauss' law for magnetic fields, for which the integral formu- 
lation is given by (1.2b). These four equations as a whole are associated with 
the name of Maxwell because he was responsible for completing them, adding 
to Ampere's original equation, V x H(f,t) = J(f,t), the displacement current 
density term or, in short, the displacement current, dD/dt, as an additional 
vector source for the field H . This term has the same dimensions as the free 
current density but its nature is different because no free charge movement is 
involved. Its inclusion in Maxwell's equations is fundamental to predict the ex- 
istence of electromagnetic waves which can propagate through empty space at 
the constant velocity of light c. The concept of displacement current is also fun- 
damental to deduce from (l.ld) the principle of charge conservation by means 
of the continuity equation 

dp 



or, in integral form, 



*-J = -f t (1-3) 



With his equations, Maxwell validated the concept of "field" previously in- 
troduced by Faraday to explain the remote interactions of charges and currents, 
and showed not only that the electric and magnetic fields are interrelated but 
also that they are in fact two aspects of a single concept, the electromagnetic 
field. 

The link between electromagnetism and mechanics is given by the empirical 
Lorenz force equation, which gives the electromagnetic force density, / (in N 
m~ 3 ), acting on a volume charge density p moving at a velocity u (in m s^ 1 ) 
in a region where an electromagnetic field exists, 

/ '= p(E + ux B) = pE + fx B (1.5) 

where J = pu is the current density in terms of the mean drift velocity of the 
particles 5 , which is independent of any random velocity due to collisions. The 



4 Given that the tesla is an excessivelly high magnitude to express the values of the magnetic 
field usually found in practice, the cgs unit (gauss, G) is often used instead, IT = 10 4 G. 

5 In general, when there is more than one type of particle the current density its defined 
as J = Pi^i where p i and u; represent the volume charge density and drift velocity of the 
charges of class i. 
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total force F exerted on a volume of charge is calculated by integrating / in this 
volume. For a single particle with charge q the Lorentz force is 

F = q(E + uxB) (1.6) 

Maxwell's equations together with Lorenz's force constitute the basic mathe- 
matical formulation of the physical laws that at a macroscopic level explain and 
predict all the electromagnetic phenomena which basically comprise the remote 
interaction of charges and currents taking place via the electric and/or magnetic 
fields that they produce. From Eq. (1.6) the work done by an electromagnetic 
field acting on a volume charge density p inside a volume dv during a time 
interval dt is 

dW = f ■ udtdv = p(E + ux B) ■ udtdv = pE ■ udtdv = E ■ Jdtdv (1.7) 

This work is transformed into heat. The corresponding power density P v ( 
WrrT 3 ) that the electromagnetic field supplies to the charge distribution is 

p _dP _ dW _ ~ j 
v dv dtdv 

This equation is known as the point form of Joule's law. 

In applications, Maxwell's equations have to be complemented by appropri- 
ate initial and boundary conditions. The initial conditions involve values or 
derivatives of the fields at t = 0, while the boundary conditions involve the 
values or derivatives of the fields on the boundary of the spatial region of inter- 
est. Usually, we consider the initial conditions as a form of boundary conditions 
and refer to the solution of Maxwell's equations, with all these conditions, as a 
boundary-value problem. 

Next, we briefly describe the physical meaning of Maxwell's equations. 



1.2.1 Physical meaning of Maxwell's equations 

Gauss' law, (1.1a) or (1.2a), is a direct mathematical consequence of Coulomb's 
law, which states that the interaction force between electric charges depends on 
the distance, r, between them, as r~ 2 . According to Gauss' law, the divergence 
of the vector field D is the volume density of free electric charges which are 
sources or sinks of the field D, i.e. the lines of D begin on positive charges 
(p > 0) and end on negative charges (p < 0). In its integral form, Gauss' 
law relates the flux of the vector D through a closed surface S (which can be 
imaginary; Fig. 1.1), to the total free charge within that surface. 

Gauss' law for magnetic fields, (1.1b) or (1.2b), states that the B field does 
not have scalar sources, i.e., it is divergenceless or solenoidal. This is because 
no free magnetic charges or monopoles have been found in nature (see Section 
2.5) which would be the magnetic analogues of electric charges for E. Hence, 
there are no sources or sinks where the field lines of B start or finish, i.e., the 
field lines of B are closed. In its integral form, this indicates that the flux of 
the B field through any closed surface S is null. 
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Figure 1.1: (a) Closed surface *S* bounding a volume V. (b) Open surface S bounded by the 
closed loop r. The direction of the surface element dS is given by the right-hand rule: the 
thumb of the right hand is pointed in the direction of dS and the fingertips give the sense of the 
line integral over the contour F. jj jAtencion: las dS debe ser ds 



Faraday's law, (1.1c) or (1.2c), establishes that a time-varying B field pro- 
duces a nonconservative electric field whose field lines are closed. In its integral 
form, Faraday's law states that the time variation of the magnetic flux (J B ■ ds) 
through any surface S bounded by an arbitrary closed loop F, (Fig. 1.1), in- 
duces an electromotive force given by the integral of the tangential component 
of the induced electric field around T. The line integration over the contour T 
must be consistent with the direction of the surface vector ds according to 
the right-hand rule. The minus sign in (1.1c) and (1.2c) represents the feature 
by which the induced electric field, when it acts on charges, would produce an 
induced current that opposes the change in the magnetic flux (Lenz's law). 



Ampere's generalized law, (l.ld) or (1.2d), constitutes another connection, 
different from Faraday's law, between E and B. It states that the vector sources 
of the magnetic field may be free currents, J, and/or displacement currents, 
dD/dt. Thus, the displacement current performs, as a vector source of H, a 
similar role to that played by dB/dt as a source of E. In its integral form the 
left-hand side of the generalized Ampere's law equation represents the integral 
of the magnetic field tangential component along an arbitrary closed loop V 
and the right-hand side is the sum of the flux, through any surface S bounded 
by a closed loop V (Fig. 1.1 ), of both currents: the free current J and the 
displacement current dD/dt. 
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1.2.2 Constitutive equations 

Maxwell's equations (1.1) can be written without using the artificial fields D 
and H, as 



V • E{r, t) 
V-B(r,t) 

V x E(r, t) 

V x B(r, t) 



^(r,i) 
so 



dB(r, t) 
dt 

H Jall{r,t) + /IqEq 



dE(f,t) 
dt 



(1.9a) 
(1.9b) 
(1.9c) 

(1.9d) 



where So = 10 _9 /(367r) (farad/meter; F to -1 ) and p, Q = 47rl0~ 7 (henry/meter; 
H m _1 ) are two constants called electric permittivity and magnetic permeability 
of free space, respectively. The subscript all indicates that all kinds of charges 
(free and bound ) must be individually included in p and J. These equations 
are, within the limits of classical electromagnetic theory, absolutely general. 
Nevertheless, in order to make it possible to study the interaction between an 
electromagnetic field and a medium and to take into account the discrete nature 
of matter, it is absolutely necessary to develop macroscopic models to extend 
equations (??) and (??) and to obtain Maxwell's macroscopic equations (1.1), 
in which only macroscopic quantities are used and in which only the densi- 
ties of free charges and currents explicitly appear as sources of the fields. To 
this end, the atomic and molecular physical properties, which fluctuate greatly 
over atomic distances, are averaged over microscopically large- volume elements, 
Aw, so that these contain a large number of molecules but at the same time 
are macroscopically small enough to represent accurate spatial dependence at a 
macroscopic scale. As a result of this average, the properties of matter related 
to atomic and molecular charges and currents are described by the macroscopic 
parameters, electric permittivity e, magnetic permeability p, and electrical con- 
ductivity a. These parameters, called constitutive parameters, are in general 
smoothed point functions. The derivation of the constitutive parameters of a 
medium from its microscopic properties is, in general, an involved process that 
may require complex models of molecules as well as quantum and statistical 
theory to describe their collective behavior. Fortunately, in most of the practi- 
cal situations, it is possible to achieve good results using simplified microscopic 
models. Appendices ?? and ?? present a brief introduction to the microscopic 
theory of electric and magnetic media, respectively. 

To define the electric permittivity and describe the behaviour of the electric 
field in the presence of matter, we must introduce a new macroscopic field 
quantity P (C m~ 2 ), called electric polarization vector, such that 



D = s E + P 



(1.10) 
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and denned as the average dipole moment per unit volume 

P= lim ^"r 1 Pn (1.11) 

Av->0 Av 

where N is the number of molecules per unit volume and the numerator is 
the vector sum of the individual dipolar moments, p n , of atoms and molecules 
contained in a macroscopically infinitessimal volume Av. For many materials, 
called linear isotropic media, P can be considered colinear and proportional to 
the electric field applied. Thus we have 

P = eoX e E (1.12) 

where the dimensionless parameter x e , called the electric susceptibility of the 
medium, describes the capability of a dielectric to be polarized. Expression 
(1.10) can be written in a more compact form as 

D = (l + Xe )e E (1.13) 

so that 

D = e e r E = sE (1.14) 

where 

Sr = l+Xe (1-15) 

and 

£ = eos r (1-16) 

are the relative permittivity and the permittivity of the medium, respectively. 

To define the magnetic permeability and describe the behaviour of the mag- 
netic field in the presence of magnetic materials, we must introduce another 
new macroscopic field quantity, called magnetization vector M (A to -1 ), such 
that 

H= — -M (1.17) 
Mo 

where M is defined, in a similar way to that of the electric polarization vector, 
as the average magnetic dipole moment per unit volume 

M = lim mn (1.18) 

Av^O Av 

where N is the number of atomic current elements per unit volume and the 
numerator is the vector sum of the individual magnetic moments, m n contained 
in a macroscopically infinitessimal volume Aw. 

In general, M is a function of the history of B or H, which is expressed 
by the hysteresis curve. Nevertheless, many magnetic media can be considered 
isotropic and linear, such that 

M = x m H (1.19) 
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where x m i s the adimensional magnetic susceptibility magnitude, being negative 
and small for diamagnets, positive and small for paramagnets, and positive and 
large for ferromagnets. Thus 

B = (1 + xJMo# = AVMo# = M# (1-20) 

where 

M r = (l + Xj (1-21) 

and 

H = fi r fj, (1.22) 

are the relative magnetic permeability and the permeability of the medium, 
respectively, which can reach very high values in magnetic materials such as 
iron and nickel. 

The concept of \i T requires a careful definition when working with magnetic 
materials with strong hysteresis, such as ferromagnetic media. The phenomenon 
of hysteresis may also occur in certain dielectric materials called ferroelectric (see 
Appendix ??). 

In a vacuum, or free space, e r = 1; fi r = 1, and therefore the fields vectors 
D and E, as well as B and H, are related by 

D = s E (1.23a) 
B = n H (1.23b) 

Very often the relation between an electric field and the conduction current 
density J c that it generates is given, at any point of the conducting material, 
by the phenomenological relation, called Ohm's law 

J c = aE (1.24) 

so that J is linearly related to E trough the proportionality factor (7 called the conductivity of the 
medium. Conductivity is measured in Siemens per meter (S m 1 = fl l m 1 ) or mhos per 

meter (mho m _1 ). Media in which (1.24) is valid are called ohmic media. A 
typical example of ohmic media are metals where (1.24) holds in a wide range 
of circumstances. However, in other materials, such as semiconductors, (1.24) 

it may not be applicable. For most metals a is a scalar with a magnitude that depends on 
the temperature and that, at room temperature, has a very high value of the order of 10 7 mho 
m _1 .Thcn very often metals are considered as perfect conductors with an infinite conductivity. 

The relations between macroscopic quantities, (1.13), (1.20) and (1.24), arc called constitutive 
relations. Depending on the characteristics of the constitutive macroscopic parameters E, fJ, and 
(J, which are associated with the microscopic response of atoms and molecules in the medium, this 
medium can classified as: 

Nonhomogcncous or homogeneous: according to whether or not the constitutive parameter of 
interest is a function of the position, £ — e{t), /I = ^{t), or (J = g{t). 

Anisotropic or isotropic: according to whether or not the response of the 
medium depends on the orientation of the field. In isotropic media all the 
magnitudes of interest are parallel, i.e., E and D; and/or E and J c ; and/or 
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B and H. In anisotropic materials the constitutive parameter of interest is a 
tensor (see Chapter ??) 

Nonlinear or linear: according to whether or not the constitutive parameters 
depend on the magnitude of the applied fields. For instance s(E), n(H) or a(E) 
en general funcion de E y Bit 

Time-invariant: if the constitutive parameters do not vary with time e ^ 
e(t), [i ^ fi(t) or a ^ <r(t) 

Dispersive: according to whether or not, for time-harmonic fields, the con- 
stitutive parameters depend on the frequency, e = e(w), fJ. = or a = <j(uj). 
The materials in which these parameters are functions of the frequency are 
called dispersive 6 . 

Magnetic medium: if [i ^ fi . Otherwise the medium is called nonmagnetic 
because its only significant reaction to the electromagnetic field is polarization. 

Fortunately, in many cases the medium in which the electromagnetic field ex- 
ists can be considered homogeneous, linear and isotropic, time-invariant, nondis- 
persive and nonmagnetic. Indeed, this assumption is not very restrictive since 
many electromagnetic phenomena can be studied using this simplification. In 
fact, even practical cases of the propagation of electromagnetic waves through 
nonlinear media (semiconductors, ferrites, nonlinear crystals, etc.) are analysed 
with linear models using the so-called small-signal approach. Most of this 
book concerns homogeneous, linear, isotropic and nonmagnetic media, except 
in Chapter ?? where anisotropic and magnetic materials (ferrites) are consid- 
ered. 

The effect of the properties of a medium on the macroscopic field can be 
emphasized by expressing E and B in Maxwell's equations (1.1a) and (l.ld) by 
(1.10) and (1.17). Thus we have 

V.£ = P^L = — ( p -X/.p) (i.25a) 

VxB = ^J a ii+^e —^^{J+ — + VxM)+e a ^ — 

(1.25b) 



Eqs (1.12), (1.19) and (1.24) are strictly valid only for nondispersive media Effectively, for 
example, because of the dependence of the electric permittivity with frequency we generally have 
P{oj) = SoXei^E^). Thus, according to the convolution theorem, for arbitrary time dependence 
this expression becomes 

P(t) = e [ Xe (t-t')E(t')dt' 

Similarly for magnetization and Ohms' law we have 

M{t) = f Xm(t-t')H{t')dt' 

J — OO 

J(t) = f a{t - t')E(t')dt' 

J —OO 

These expressions indicate that, as for any physical system, the response of the medium to an 
applied field is not instantaneous. 



12 CHAPTER 1. ELECTROMAGNETIC FIELD FUNDAMENTALS 



In (1.25a) we have explicitly as scalar sources of E both the free charge p 
and the polarization or bounded density of charge, —V • P. Then in (1.9a) we 
have 

P all =p-V-P (1.26) 

Similarly, in (1.25b), we have, explicitly as vector sources of B, besides the free current density 
J ( which includes the conduction current density J c — <jE), the polarization Current 

dP/dt (which results from the motion of the bounded charges in dielectrics), the 
displacement current in the vacuum, SodE/dt and the magnetization current, 
V x M (which takes place when a non-uniformly magnetized medium exists). 
Then in (1.9d) we have 

- dP 

Jaii(r 7 t) = J+ — + VxM (1.27) 

In the following we will assume that there is no magnetization current. 



1.2.3 Boundary conditions 

As is evident from (l.la)-(l.ld) and (1.13), (1.20), (1.24), in general the fields 
E, B, D and H are discontinuous at points where e, p and a also are. Hence 
the field vectors will be discontinuous at a boundary between two media with 
different constitutive parameters. 

The integral form of Maxwell's equations can be used to determine the 
relations, called boundary conditions, of the normal and tangential components 
of the fields at the interface between two regions with different constitutive 
parameters e, p and er where surface density of sources may exist along the 
boundary. 

The boundary condition for D can be calculated using a very thin, small pill- 
box that crosses the interface of the two media, as shown in Fig. 1.2. Applying 
the divergence theorem 7 to (1.1a) we have 



<j)D.ds= J D\.ds+ J D.ds+ J D 2 .ds = j pdv (1.28) 

Base 1 Curved surface Base 2 

where D\ denotes the value of D in medium 1, and D 2 the value in medium 
2. Since both bases of the pillbox can be made as small as we like, the total 
outward flux of D over them is (D n \ — D n2 )ds — {D\ —D^)- hds, where these D n 
are the normal components of D, ds is the area of each base, and h is the unit 
normal drawn from medium 2 to medium 1. At the limit, by taking a shallow 
enough pillbox, we can disregard the flux over the curved surface, whereupon 
the sources of D reduce to the density of surface free charge p s on the interface 

n ■ (Di - D 2 ) = p s (1.29) 



7 E1 tcorcma dc la divergencia requiere que las propiedad.es del medio van'en dc forma contmua, 
pero puede suponerse una transicion rapida pero contfnua del medio 1 al 2 
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Pilbox Infinitesimal loop 



Medium 1 
Medium 2 



4dh 



dl 



Figure 1.2: Derivation of boundary conditions at the interface of two media. 
Pintar solo la h hacia arriba y las ds una Hcia arriba y la de abajo hacia abajo 
Cuidado pilbox es con dos 1 



Hence the normal component of D changes discontinously across the interface by 
an amount equal to the free charge surface density p s on the surface boundary. 

Similarly the boundary condition for B can be established using the Gauss' 
law for magnetic fields (1.1b). Since the magnetic field is solenoidal, it follows 
that the normal components of B are continuous across the interface between 
two media 

n ■ {Bi - B 2 ) = (1.30) 

The behavior of the tangential components of E can be determined using 
a infinitesimal rectangular loop at the interface which has sides of lengh dh, 
normal to the interface, and sides of lengh dl parallel to it (Fig. 1.2). From 
the integral form of the Faraday's law, (1.2c) and defining i as the unit tangent 
vector parallel to the direction of integration on the upper side of the loop, we 
have 

(Ei -i — Ei -t)dl+ contributions of sides dh 

In the limit, as dh — > 0, the area ds — dldh bounded by the loop approaches 
zero and, since B is finite, the flux of B vanishes. Hence(£i — E 2 ) ■ i = 
and we conclude that the tangential components of E are continuous across the 
interface between two media. In terms of the normal n to the boundary, this 
can be written as 

h x (Ex - E 2 ) = (1.32) 

Analogously, using the same infinitesimal rectangular loop, it can be deduced 
from the generalized Ampere's law, (1.2d), that 

(Hi -i — H 2 ■ t)dl + contributions of sides dh 

= -(§ + j)-ds (1.33) 
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where, since D is finite, its flux vanishes. Nevertheless, the flux of the surface 
current can have a non-zero value when the integration loop is reduced to zero, 
if the conductivity a of the medium 2, and consequently J s , is infinite. This 
requires the surface to be a perfect conductor. Thus 

h x (i?! - H 2 ) = J s (1.34) 

the tangential component of H is discontinuous by the amount of surface current 
density J s . For finite conductivity, the tangential magnetic field is continuous 
across the boundary. 

A summary of the boundary conditions, given in (1.35), are particularized 
in (1.36) for the case when the medium 2 is a perfect conductor {a 2 — * oo). 



(1.35a) 
(1.35b) 
(1.35c) 
(1.35d) 

Boundary conditions when the medium 2 is a perfect conductor {a 2 — * oo) 

h x Ex = (1.36a) 

h x Hi = J s (1.36b) 

n-Di = p s (1.36c) 

h-Bi = (1.36d) 



n x {Ei 


-E 2 ) 




= 


n x {Hi 


-H 2 ) 


= Js 


n-(Di 


-D 2 ) 


= Ps 


h-{Bi 


-B 2 ) 


= 



1.3 The conservation of energy. Poynting's the- 
orem 

Poynting's theorem represents the electromagnetic energy-conservation law. To 
derive the theorem, let us calculate the divergence of the vector field Ex H in a 
homogeneous, linear and isotropic finite region V bounded by a closed surface S. 
If we assume that V contains power sources (generators) generating currents 
J, then, from Maxwell's equations (1.1c) and (l.ld), we get 

r9R f)f) 

y-{ExH)=H-VxE-E-VxH = —H ■ — E ■ — E ■ {aE + J) (1.37) 

at at 

where J represents the source current density distribution which is the primary 
origin of the electromagnetic fields 8 , while the induced conduction current den- 
sity is written as J c = aE (1-24). 



The source current may be maintained by external power sources or generators (this current is 
often called driven or impressed current). 
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As the medium is assumed to be linear, the derivates with respect to time 
can be written as 

=. dD - dE d (\ 2 \ d (\ - -\ 

E ■ — = sE- — = — [ -sE 2 = — \-E ■ D) 1.38a 
dt dt dt\2 J dt\2 J v ' 

By introducing the equalities (1.38a) and (1.38b) into (1.37), integrating over 
the volume V, applying the divergence theorem, and then rearranging terms, 
we have 

/ J-Edv = -^- f \{E-D + B-H)dv- f uE 2 dv - (f (E x H) ■ ds (1.39) 

JV Ot Jy 2 Jy Jg 

To interpret this result we accept that 

U ev = \d-E (1.40) 

and 

U mv = l -B-H (1.41) 

represent, as a generalization of their expression for static fields, the instanta- 
neous electric energy density, U ev , and magnetic energy density, U mv , stored in 
the respective fields. Thus according to (1.8) the left side of (1.39) represents 
the total electromagnetic power supplied by all the sources within the volume 
V. Regarding the right side of (1.39), the first term represents the change rate 
of the stored electromagnetic energy within the volume; the second term repre- 
sents the dissipation rate of electromagnetic energy within the volume; and the 
third term represents the flow of electromagnetic energy per second (power) 
through the surface S that bounds volume V. Defining Poynting's vector V as 

V = ExH (W/m 2 ) (1.42) 

we can write 

j>{ExH)-ds = j> V -ds (1.43) 

This equation represents the total flow of power passing through the closed sur- 
face S and, consequently, we conclude that V = E x H represents the power 
passing through a unit area perpendicular to the direction of V . This conclu- 
sion may seem questionable because it could be argued that any vector with an 
integral of zero over the closed surface S could be added to V without affecting 
the total flow. Nevertheless, this is a natural interpretation that does not con- 
tradict any experience. Only when we try to particularize (1.39) to steady fields 
do we find ambiguous results, because, in static, the location of the electric and 
magnetic energy has no physical significance. 
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Note that Eq. (1.39) was deduced by assuming a linear medium and that the 
losses occur only through conduction currents. Otherwise the equation should 
be modified to include other kinds of losses such as those due to hysteresis or 
possible transformations of the electromagnetic energy into mechanical energy, 
etc. When there are no sources within V, (1.39) represents an energy balance 
of that flowing through S versus that stored and dissipated in V. 



1.4 Momentum of the electromagnetic field 

As we have seen in the previous section, when we apply the law of conservation 
of electromagnetic energy to a finite volume V bounded by a surface S, it is 
necessary to include a term that, by means of the Poynting vector V, takes 
into account the flow of power through S. We shall now see that when an 
electromagnetic field interacts with the charges and currents in V, it is also 
necessary to consider a momentum associated with the electromagnetic field in 
order to guarantee the conservation of momentum. To calculate this momentum, 
we will begin by expressing, only in terms of the fields, the Lorentz force density, 
(1.5), exerted by the electromagnetic field on the distribution of charges and 
current, which we assume to be in free space. For this purpose, let us consider 
Maxwell's equations (1.1a) and (l.ld) to express p and J as 

p = VD (1.44) 
-. dD 

J = ^xH- — (1-45) 



so that 



f = pE + JxB=(s7-D^E-Bx(VxH) + Bx^ (1.46) 



which, taking into account that 



-^0 x B) - D x (V x E) (1.47) 



becomes 



/= (V • D)E - B x (V x H) - ^(D x B) - D x (V x E) (1.48) 

By adding the term H (V • B) = to this equality to make the final expres- 
sions symmetrical, and by reordering, we can write the Lorentz force density 
as 

/ = E(\/ ■ D) - D x (V x E) + HV ■ B - B x (V x H) - ^-{D x B) (1.49) 
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The component a of Lorentz force density can be written, taking into account 
the definition of the Poynting vector P, as 



fa 



- A 

dp 



EjiE a - -5p a E 2 



d 



'P L 



Hp Ha - -j&liuH 2 



l__d 



P a (1.50) 



where Sp a is the Kronecker delta (Spa = 1 if (5 — a and zero if /? ^ a) and the 
indices a, (3 = 1,2,3 correspond to the coordinates x, y, z, respectively, and we 
have made use of the Einstein's summation convention (i.e., the repetition of 
an index automatically implies a summation over it). To obtain (1.50) we have 
made use of the following equalities 



E a \7 ■ D- Dx (V x E) 
B a V ■ B - B x (V x H) 
Dx B 



d 



Pa 
n 2 



d 



E$E a - -5p a E 
HpH a - -5p a H 2 



(1.51) 



The first two summands in (1.50) constitute the a component of the diver- 
gence of a tensor quantity, T em , such that 



(V • T em ) c 



01 Pa 



(1.52) 



where T em is a symmetric tensor, known as the Maxwell stress tensor, defined 

by 



rpem 

1 Pa — fc o 



E^Ea^-SpaE 2 



HpH a - -Sp a H 



Therefore, from (1.50) and (1.52), we have 



\_&P_ 



with 



V-T e 



d d d 
dx' <9y ' dz 



rpem rpem 

xx xy 

rpem rpem 

± yx ± yy 

rpem rpem 

zx zy 



rpem 

xz 
rpem 

yz 

rpem 



The components of the electromagnetic tensor can be written as 



(1.53) 



(1.54) 



(1.55) 



rpem rpe , rpm r\ rp 

1 Rn — 1 Ra + 1 Ra — U P^a — 



^E^D^ + B»H a — 



(1.56) 



where and T| Q represent, respectively, the electric and magnetic tensors 
defined by 



T5 



DpE a — -5j3 a E 7 D 7 



Bf}H a 



3 a H 1 B 1 



(1.57) 
(1.58) 
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Integrating (1.50) over the volume V the total electromagnetic force F ex- 
erted on the volume is 

F = [ fdv = [ ( P E + JxB)dv= f f a da-\%- f Vdv (1.59) 
Jv Jv Js c vt Jv 

where f s is the force per unit of area on S 

f s =T em -n (1.60) 

and we have applied the theorem of divergence to the tensor T em i.e. 

/ V • T em dv = [ T em -ds= [ T em -nds= [ f s ds (1.61) 
Jv Js Js Js 

Thus 

Jj a ds = F + ^l t J v Vdv (1.62) 

Note that the term 




(1.63) 



is not null even in the absence of charges and currents. Since the only elec- 
tromagnetic force possible due to the interaction of the field with charges and 
currents is F, the term (1.63) must represent another physical quantity with 
the same dimensions as a force, i.e., the rate of momentum transmitted by 
the electromagnetic field to the volume V. This is equivalent to associating a 
momentum density g with the electromagnetic field, given by 1/c 2 times the 
Poynting vector, 

which propagates in the same direction as the flow of energy. Thus, Eq. 1.62 
represents the formulation for the momentum conservation in the presence of 
electromagnetic fields. 

The momentum of an electromagnetic field, which can be determined ex- 
perimentally, is inappreciable under normal conditions and its value is often 
below the limits of the measurement error. However, in the domain of atomic 
phenomena, the momentum of an electromagnetic field can be comparable to 
that of particles, and plays a crucial role in all the processes of interaction with 
matter. The transfer of momentum to a system of charges and currents implies 
a reduction in the field momentum, and the loss of momentum by the system, 
for example by radiation, leaves to an increase in the momentum of the field. 



1.5 Time- harmonic electromagnetic fields 

A particular case of great interest is one in which the sources vary sinusoidally 
in time. In linear media the time-harmonic dependence of the sources gives rise 
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to fields which, once having reached the steady state, also vary sinusoidally in 
time. However, time-harmonic analysis is important not only because many 
electromagnetic systems operate with signals that are practically harmonic, but 
also because arbitrary periodic time functions can be expanded into Fourier 
series of harmonic sinusoidal components while transient nonperiodic functions 
can be expressed as Fourier integrals. Thus, since the Maxwell's equations are 
linear differential equations, the total fields can be synthesized from its Fourier 
components. 

Analytically, the time-harmonic variation is expressed using the complex 
exponential notation based on Euler's formula, where it is understood that the 
physical fields are obtained by taking the real part, whereas their imaginary 
part is discarded. For example, an electric field with time-harmonic dependence 
given by cos(w£ + ip), where uj is the angular frequency, is expressed as 

E = Re{Ee jut } = ^(Ee jut + (Ee jut )*) = E cos(ojt + ip) (1.65) 

where E is the complex phasor, 

E = E e jv (1.66) 

of amplitude Eq and phase p, which will in general be a function of the angular 
frequency and coordinates. The asterisk * indicates the complex conjugate, 
and Re{} represents the real part of what is in curly brackets. 

Throughout the book, we will represent both complex phasor magnitudes 
(either scalar or vector) by symbols in bold, e.g. E — E{r,w), and p = 
p(f,oj). In this way, time-dependent (real) quantities, which are represented by 
mathematical symbols not in bold, such as E = E(r,t), and p — p(r,t), can be 
distinguished from complex phasors which do not depend on time. In general, 
as indicated, these complex phasors may depend on the angular frequency. The 
real time-dependent quantity associated with a complex phasor is calculated, as 
in (1.65), by multiplying it by e- 7 "* and taking the real part. 

1.5.1 Maxwell's equations for time-harmonic fields 

Assuming e 3 "* time dependence, we can get the phasor form or time-harmonic 
form of Maxwell's equations simply by changing the operator d/dt to the factor 
juj in (l.la)-(1.2d) and eliminating the factor e J "*. Maxwell's equations in 
differential and integral forms for time-harmonic fields are given below. 

Differential form of Maxwell's equations for time-harmonic fields 



V D = p (Gauss' law) (1.67a) 

V • B = (Gauss' law for magnetic fields) (1.67b) 
VxE = -ju B (Faraday's law) (1.67c) 
VxH = J+ jujD (Generalized Ampere's law) (1.67d) 



Integral form of Maxwell's equations for time harmonic fields 
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J D ■ ds = Q T (Gauss' law) (1.68a) 

<j> B ■ ds = (Gauss' law for magnetic fields) (1.68b) 

<j>E-dl = —ju J B ■ ds (Faraday's law) (1.68c) 

<j> H ■ dl = J (J + jtoD) ■ ds (Generalized Ampere's law) (1.68d) 
For time-harmonic fields, expressions (1.25a) and (1.25b) become 

V • E = P^L = — (p-y.p) (i.69a) 

VxB = jweoMo ^ + A*o Jail = jujeqUqE + /i (J + jujP + V x M) 

(1.69b) 

1.5.2 Complex dielectric constant. 

Over certain frequency ranges, due to the atomic and molecular processes in- 
volved in the macroscopic response of a medium to an electromagnetic field, 
there appear relatively strong damping forces that give rise to a delay between 
the polarization vector P and E (a phase shift between P and E), and con- 
sequently between E and D, and to a loss of electromagnetic energy as heat in 
overcoming the damping forces (see Appendix ??). At the macroscopic level this effect 
is analytically expressed by means of a complex permittivity. £ c as 

D = e c E (1.70) 

with 

e c = e' - js" = e e cr (1-71) 

where 6 cr 

£cr = l+Xce=4-i4' (1-72) 
is the relative complex permittivity and \ ce = Xcer — JXcer ^ s t ^ lc com pl° x electric susceptibility. 

In general both e' and e" present a strong frequency dependence and they 
are closely related to one another by the Kramer-Kronig relations as is shown 
in Appendix ??, where the dependence with the frequency of the dielectric 
constant is studied. 

Similar processes occur in magnetic and conducting media, and, within a 
given frequency range, there may be a phase shift between E and J c or between 
B and H which, at the macroscopic level, is reflected in the corresponding 
complex constitutive parameters er c = a' — ja" and /j, c = fjf — jfi". 

For a medium with complex permittivity, the complex phasor form of the 
displacement current is 

jujD = jujs c E = ue"E + jcue'E (1.73a) 
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Figure 1.3: Induced current density in the complex plane. 

while the sum, of the displacement and conduction current, called total induced 
current, J$, is 

J* = aE + ju>e c E = (a + cue")E + jcus'E = J d + J r (1.74) 

where Jd, called the dissipative current, 

J d = (a + Loe")E (1.75) 

in phase with the electric field, is the real part of the induced current J i (Fig. 
1.3) while J r , called the reactive current, 

J r = jcus'E (1.76) 

is the imaginary part of the induced current which is in phase quadrature with 
the electric field. The dissipative current can be expressed in a more compact 
form as 

Jd = o e E (1.77) 
where er e is the effective or equivalent conductivity 

a e = a + ue" (1.78) 

which includes the ohmic losses due to a and the damping losses due to we". 
Thus the induced current, (1-74), can be written as 

J* = a e E + jcus'E = a ec E (1.79) 

where er ec is the complex effective conductivity defined as 

<r ec = <j e +juje' (1.80) 

Thus a medium with conductivity a ec and null permittivity is formally equiva- 
lent to one with conductivity and permittivity, a and e c , respectively. 
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On the other hand, the phase angle 8d between the induced and reactive 
currents, (Fig. 1.3), is called the loss or dissipative angle, and its tangent (i.e., 
the ratio of the dissipative and reactive currents) is called the loss tangent 

tan<5 d = ^- (1.81) 
we' 

and the induced current, (1.79), can be written in terms of the loss tangent as 

J* = a ec E — jus' (I - j^j)E = jue' (1 - j temS d ) E = jcoe ec E (1.82) 

where e ec is defined as the effective complex permittivity 

e ec = e' (1 - jt&nSd) = e e er (1.83) 

and 

e er = (1 -jtanS d )e' r (1.84) 

denotes the effective relative permittivity. Thus, according to (1.79) and (1.82), a medium 
can be formally considered alternatively cither as a medium of permittivity &' and effective con- 
ductivity <7 e , or as a dielectric medium of effective permittivity £ ec or as a conducting medium of 
effective conductivity (7 ec . In summary, this possibilities arc 



Permittivity Conductivity 
Original medium e c = e' — je" a 

Equivalent medium 1 e' a e = a + cue" 

Equivalent medium 2 e ec = e' — j(s" + ^) 
Equivalent medium 3 cr ec = a + ue" + jue' 

(1.85) 

The loss tangent is equal to the inverse of the quality factor Q of the dielectric which is a 
dinicnsionlcss quantity defined as 



Maximun energy stored per unit volume W v 

Q = Ul = U! 

Time average power lost per unit volume P dv 



(1.86) 



The average power dissipated per cycle and unit volume, P' dv , due both to the 
Joule effect and to that of dielectric polarization, is given, according to (1.8), 

by 

1 f T ^ 1 f T 

P' dv = — I E ■J i dt=— I E coswi • (cr e Eo cosut + uie'Eq s'mujt)dt 
1 Jo 1 Jo 

1 1 G E^ 

= - / a e E 2 cos 2 ujtdt =- E ■ J d dt = 

(1.87) 
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where T = 2n/ui is the period of the signal. Note that only the dissipative part 
of Ji contributes to the average power. Of this power, the part corresponding 
to polarization losses is 

^ f T us"E 2 cos 2 ujtdt = (1.88) 
*■ Jo 2 

The maximum electric field energy stored per unit of volume is 

W v = l -e'El (1.89) 
Thus, dividing (1.89) by (1.87), we have 

Q = — = r—T- (1-90) 

a e tandrf 

Although both dimensionless quantities, Q and tan^, can be used to define 
the characteristics of a dielectric, we will use the loss tangent throughout this 
book. 

Depending on whether the reactive or the dissipative current is predominant 
at the operating frequency, a medium is classified as a weakly lossy or a strongly 
lossy medium respectively. Thus for weakly lossy media, usually called good 
dielectrics or insulators, we have, us' » a e , so that 

tanJ d =-^-«l (1.91) 
we' 

Or, if a = 0, 

e" 

t&n5 d = — «l (1.92) 

s 

If <7 e = (i.e. tan<5d = 0), the medium is termed a perfect or ideal dielectric, 
in which case the reactive current coincides with the displacement current, and 
the dielectric is characterized by a real permittivity e. 

If the medium is strongly lossy we have cue' « a e , so that 

tan5 d =^»l (1.93) 

which for good conductors where e" = 0; e' = e simplifies to 

tan5 d = — >> 1 (1.94) 

cue 

being practically s = eq. If a = oo (i.e. tan<5 d = oo) the medium is termed a 
perfect conductor. 

For a homogeneous conducting medium where e' and a e do not depend on 
the position, Gauss' law (1.1a) and the continuity equation (1.3) can be writen 
as 

V • E = p/e' (1.95) 
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and 

<7 e V-E = -^ (1.96) 

respectively. Hence we have 

so that the expression for the decay of a charge distribution in a conductor is 
given by 

p = p e-(^'> (1.98) 
where p is the charge density at time t = 0. The characteristic time 

r=— (1.99) 

required for the charge at any point to decay to 1 /e of its original value is called 
the relaxation time. 

For most metals r = lCT 14 s, signifying that in good conductors the charge 
distribution decays exponentially so quickly that it may be assumed that p = 
at any time. In terms of the relaxation time, the loss tangent can be written as 

tan 5 d = — = (tuj)- 1 (1.100) 
ew 

Thus the classification of a medium as a good or poor conductor depends on 
whether the relaxation time is short or long compared with the period of the 
signal. 

1.5.3 Boundary conditions for harmonic signals 

For harmonic signals the boundary conditions of the normal and tangential 
components of the fields at the interface between two regions with different 
constitutive parameters s, p, and a, (1.35a)-(1.36d), become 
General boundary conditions 

h x (Ei - E 2 ) = (1.101a) 

nx{H 1 -H 2 ) = f s (1.101b) 

n-(Di-D 2 ) = Ps (1.101c) 

n-(B 1 -B 2 ) = (l.lOld) 

Boundary conditions when the medium 2 is a perfect conductor (er 2 — * oo) 

n x Ei = (1.102a) 

n x Hi = J s (1.102b) 

n-Dx = p s (1.102c) 

n-Bi = (1.102d) 
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1.5.4 Complex Poynting vector 

In formulating the conservation- energy equation for time-harmonic fields, it is 
convenient to find, first, the time-average Poynting vector over a period, i.e. the 
time-average power passing through a unit area perpendicular to the direction 
of V. From (1.65) we have 



E = RejEV"*} = X - (Ee jut + (eV w *)*) (1.103a) 

H = Re{He jw *} = X - (h^' + (He iut ) j (1.103b) 
Thus, the instantaneous Poynting vector (1.42) can be written as 
V = ExH = Re{Ee JaJ '} x Re{He jujt } 

= i Re{E x H* + E x He 2jujt } (1.104) 

where we have made use of the general relation for any two complex vectors A 
and B 

Re{A} x Re{B} = i(A + A) x i(B + B) 

= -(Ax B* + A x B) + -(Ax B + A x B*) 

4 4 

= ^Ax B* + (Ax B*y^ +^Ax B + (Ax B*y^ 

= i Re{A x B* + A x B} (1.105) 

The time-average value of the instantaneous Poynting vector can be calcu- 
lated integrating (1.104) over a period , i.e., 

-| f-T i r-T 



1 f 1 f 

- pdt=— Re{E xH* +Ex He 2jujt }dt 
T Jo 2T J 



T 

= ^Re{ExH*} = ^Re{P c } (1.106) 

since the time average of E x He 2)u ' vanishes. The magnitude 

V c = ExH* (1.107) 

is termed the complex Poynting vector. Thus the time-average of the Poynting 
vector is equal to one-half the real part of the complex Poynting vector 

For a more complete view of the meaning of the complex Poynting vector, 
let us again formulate Poynting's theorem particularized for sources with time- 
harmonic dependence. From Faraday's law, (1.67c), and from Ampere's general 
law, (1.67d), in its conjugate complex form, we have 

VxE = -juj^H (1.108a) 
VxH = -jusE* +J*+ aE* (1.108b) 
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where J represents the complex conjugate of the current supplied by the 
sources. Performing a scalar multiplication of Eq. 1.108a by H and of Eq. 
1.108b by E, and subtracting the results, we get 

V^Exffj = jT-VxE-E-VxH* 

= -juj (m#o - eE l) - E ■ [J* + aE*) (1.109) 

where it has been taken into account that H ■ H = Hq and E ■ E = Eq, 
with Hq and Eq being the amplitude of the two harmonic fields. After dividing 
(1.109) by 2 we get 

V • ( U x ii*) = -2 ju (^f-S^-^-'-r-E (1.110) 



2 2 

The terms hHq/A and sEq/A represent, respectively, the mean density of the 
magnetic and electric energy while (jEq/2 is the the mean power transformed 
into heat 9 within V, since the mean value of the square of a sine or cosine 
function is 1/2. 

By multiplying Equation (1.110) by the volume element dv, integrating over 
an arbitrary volume V and applying the divergence theorem, we obtain the 
complex version of the Poynting theorem 



X 



5 ('•■*)* - -X^-^X^-t)* 

-J^(SxH*yds (1.111) 

which is the expression corresponding to (1.39) in complex notation and where 
the first member represents the power supplied by external sources. By sepa- 
rating the real and imaginary parts, we obtain the following two equalities 

[Rel(f*-E)dv = -[ °^-dv- [ Re^(E x H*) ■ ds (1.112a) 
Jv 2 J v 2 J s 2 

f lm\{J * -E)dv = -2w / (^-e^f]dv-[ljn]-{ExH)-ds 
Jv 2 Jv \ 4 4 J J s 2 

(1.112b) 

The first member of (1.112a) 

P a = I Rehf*-E)dv (1.113) 
Jv z 

represents the active mean power supplied by all the sources within V. On the 
right-hand side of (1.112a) the first integral, as commented above, gives the 



9 Expression (1.112b) can be easily extended to the case of lossy dielectric just substituting a by 
the equivalent conductivity cr e denned in (1.78)and e by e defined in (1.71). 
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power transformed into heat within V, while the surface integral represents the 
mean flow of power through the surface S. 

Regarding to expression (1.112b), the first member 



is called the reactive power of the sources. On the right-hand side the first 
summand is 2lu times the difference of the average energies stored in the 
electric and magnetic fields, while the second represents the flow of reactive 
power that is exchanged with the external medium through S. If the surface 
integral in (1.112a) is non-zero, the external region is said to be an active charge 
for the sources within V. Similarly, if the surface integral of Eq. (1.112b) is non- 
zero, the external region is said to be a reactive charge for the sources within V. 
In general, both of these surface integrals are non-zero and the external region 
becomes both an active and a reactive charge for the sources. 

1.6 On the solution of Maxwell's equations 

Despite their apparent simplicity, Maxwell's equations are in general not easy 
to solve. In fact, even in the most favorable situation of homogeneous, linear 
and isotropic media, there are not many problems of interest that can be analyt- 
ically solved except for those presenting a high degree of geometrical symmetry. 
Moreover, the frequency range of scientific and technological interest can vary 
by many orders of magnitude, expanding from frequency values of zero (or very 
low) to roughly 10 14 Hertz. The behavior and values of the constitutive para- 
meters can change very significantly in this frequency, range. Conductivity, for 
example, can vary from to 10 7 S m _1 . It is even possible to build artificial 
materials, called metamaterials, which present electromagnetic properties that 
are not found in nature. Examples of such as metamaterials are those char- 
acterized with both negative permittivity (e < 0) and negative permeability 
(fi < 0). These media are called DNG (double- negative) metamaterials and, 
owing to their unusual electromagnetic properties, they present many potential 
technological applications. 

Another important factor to study the interaction of an electromagnetic field 
with an object is the electrical size of the body, i.e., the relationship between 
the wavelength and the body size, which can also vary by several orders of mag- 
nitude. All these circumstances make it in general necessary to use analytical, 
semi-analytical or numerical methods appropriate to each situation. In partic- 
ular, numerical methods are fundamental for simulating and solving complex 
problems that do not admit analytical solutions. Today numerical methods 
make up the so-called computational electromagnetics, which together, with ex- 
perimental and theoretical or analytical electromagnetics, constitute the three 
pillars supporting research in Electromagnetics. Of course, both the develop- 
ment of analitycal, numerical or experimental tools, as well as the interpretation 
of the results, require theoretical knowledge of electromagnetic phenomena 




(1.114) 
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Chapter 2 



Fields created by a source 
distribution: retarded 
potentials 

In this chapter, we introduce the scalar electric and magnetic vector potentials 
as magnitudes that facilitate the calculation of the fields created by a bounded- 
source distribution, paying special attention to the radiation field. Finally, we 
extend Maxwell's equations, in order to make them symmetric, by introducing 
the concept of magnetic charges and currents. 



2.1 Electromagnetic potentials 

A basic problem in electromagnetism is that of finding the fields created for a 
time- varying source distribution of finite size, which we assume to be in a non- 
magnetic, lossless, homogeneous, time-invariant, linear and isotropic medium. 
Figure 2.1 represents such a distribution, where, as usual, the coordinates asso- 
ciated with source points, J = J(f,t'), p = p(r f ,t / ), are designated by primes, 
while those associated with field points or observation points P(f, t) are without 
primes. In the following, we will assume the medium surrounding the source 
distribution to be free space, i.e. fj, = /i , e = Eo, although of course all the re- 
sulting formulas remain valid for media of constant permittivity and permeabil- 
ity, provided that e is replaced by e r e and [i by \i r [i {) . While the expressions 
for the fields can be derived directly from their sources, the task can often be 
facilitated by calculating first two auxiliary functions, the scalar electric poten- 
tial <E> = $(r, t) and the magnetic vector potential A = A(f,t) (Fig. 2.2). Once 
the potentials are obtained, it is a simple matter to calculate the fields from 
them. In this section, we formulate the general expressions for these potentials. 
Since, according to (1.1b), the divergence of the magnetic field B is always 
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Figure 2.1; Time-varying source bounded distribution V' of maximun dimension I. The coor- 
dinates associated with source points of currents and charges J — J(T ^if\ and p — £ ), 
respectivelly are designated by primes, while the associate with field points, P{v 7 t^, are without 
primes. 



p(r\t)J(r\t) 


■> 


E;B 










i 












> 







Figure 2.2: Quitar o poner argumentos pero unificar 



2.1. ELECTROMAGNETIC POTENTIALS 



31 



zero, we can express it as the curl of an electromagnetic vector potential A as 

B = Vxl (2.1) 
Inserting this expression into (1.1c) we get 

Vx(i+^j=0 (2.2) 

Since any vector with a zero curl can be expressed as the gradient of a scalar 
function $, called the scalar potential, we can write 

E+^A = -V<5> (2.3) 

or 

r)A 

E = -V$ - -g (2.4) 

where dA/dt is the nonconservative part of the electric field with a non- vanishing 
curl. When the vector potential A is independent of time, expression (2.4) 
reduces to the familiar E(r) = — V$(r). 

According to the relations (2.1) and, (2.4) the fields B and E are completely 
determined by the vector and scalar potentials A and <£. However, the fields 
do not uniquely determine the potentials. For instance, it is clear that the 
transformation 

A = A' + W (2.5) 

where <]/ = ^ (r , t) is any arbitrary, single- valued, continuously differentiable, 
scalar function of position and time that vanishes at infinity, leaves B unchanged 

B = Vxl=Vxl + VxV$ = VxI' (2.6) 

Inserting (2.5) into (2.4), it follows that 

so that the value of E, obtained from A', also remains unchanged provided that 
$ is replaced by the scalar potential 

$' = $+^- (2.8) 

Thus different sets of potentials A and $ give rise to the same set of fields 1 B 
and E. The joint transformation (2.5) and (2.8) leaves the electromagnetic field 



lr The liberty to select the value of A is understandable taking into account that by (2.1) the 
magnetic field fixes only V X A. However, Helmholtz's theorem posits that, to determine the 
(spatial) behavior of A completely, V.A (which is still undetermined) must also be specified. 
Thus, we can choose it in any way we consider suitable for facilitating the calculation of the 
electromagnetic fields. 
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invariant. The different forms of choosing the potentials A and <3> leaving the 
fields unchanged are called gauge transformations, and the function \t is called 
the gauge function. The degree of freedom provided by the gauge transforma- 
tions facilitates the calculation of the potentials and hence of the fields because, 
once the potentials are known, the fields are easily derived by differentiation 
from (2.1) and (2.4). An example of gauge transformation is the Lorenz gauge, 
also called the Lorenz condition. 



2.1.1 Lorenz gauge 



Inserting (2.1) and (2.4) into the generalized Ampere's law, (l.ld), and Gauss' 
law, (1.1a), using (??) and rearranging terms, we get two, coupled, second-order 
partial-differential equations 



d2 A ^ 2 t 
War-* 2 * 



= A*o J - v 



P_ 



a_ 
dt 



V • A + ^ e 

V • A + fj, e 



<9$ 
~dt 
9$' 
~dt 



(2.9a) 
(2.9b) 



These equations could be considerably simplified if we could force (without 
changing the fields) the potentials to satisfy the auxiliary relation 



(2.10) 



called the Lorenz gauge (or Lorenz condition) 2 . Fortunately, as we will show 
below, we can always take advantage of the freedom in choosing the potentials 
so that they fulfil the Lorenz condition and consequently simplify Eqs (2.9) to 
the inhomogeneous Helmholtz wave equations 



V-A + Mo£o ^ = 



/j, J (2.11a) 
(2.11b) 

£o 

The advantage of having applied the Lorenz condition is that the equations 
(2.11) for the potentials are uncoupled and each one depends on only one type 
of source. This makes it easier to calculate the potentials than the fields (see 
Section ??). 

It remains to be shown that it is always possible to force the potentials to 
satisfy the Lorenz condition (2.10). To this end, let us consider two potentials, 
A' and which fulfil Equations (??) and (??) and check whether it is possible 



92 A V72 X 

M ° £o w~ = 

d2 $ v,2^ 



A very interesting property of the Lorenz condition is that, as shown in (??), it is covariant. i.e., 
if it holds in one particular inertial frame then it automatically holds in all other inertial frames. 



2.1. ELECTROMAGNETIC POTENTIALS 



33 



to select them so that they satisfy equations (2.11). By inserting (2.5) and (2.8) 
into (2.9) and by rearranging, we get 



d 2 A' - /„ T , ~ d& a 2 *\ 

V 2 i4'-A*oeo-^- = -Mo^ + Vf V ■ A' + + » £ — - ^0-^ ) 



(2.12a) 



(2.12b) 

and, given that the scalar function if? is arbitrary, we can choose it as the solution 
to the differential equation 

d 2 ^ - <9$' 

V 2 V - ^ s — = -V ■ A' - ti a s — (2.13) 

Thus (2.12a) and (2.12b) become (2.11a) and (2.11b), respectively, meaning 
that A' and fulfil the Lorenz condition. 

Expressions (2.11a) and (2.11b) are the inhomogeneous wave equations for 
the potentials, and their solutions, which are provided in the next section, rep- 
resent waves propagating at the velocity c = l/ A /Ai £ o — 3 x 10 8 m/s of light in 
free space. They take the form 

V ^-^^ = nA=-ft>J (2.14a) 
V 2 *-l|| = □* = -£ (2.14b) 
where the symbol □ represents the D'Alembertian operator defined by 

° <2 - 15 > 

The Lorenz gauge (2.10) for harmonic fields simplifies to 

V-A+^* = (2.16) 



such that 

while (2.14a) and (2.14b) simplify to 

j 2 
c 



(2.17) 



V 2 A+^-A = - Mo J (2.18a) 

2 

V 2 * + ^* = — ^ (2.18b) 

c 2 e 
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In addition to Lorenz's gauge, other gauge conditions may sometimes be 
useful. For instance, in quantum field theory, where the potentials are used to 
describe the interaction of the charges with the electromagnetic held instead of 
being used to calculate the helds, it is useful to use Coulomb's gauge, in which 
V • A = 0. By taking the divergence of (2.4) and the curl of (2.1), and taking 
into account the generalized Ampere's law and Gauss' law , we can easily see 
that with Coulomb's gauge the expressions for the potential wave equations are 

v 2 $ = -L 

£o 

2 r 1 d 2 A 1 a$ T 

As can be seen from (2.19), in Coulomb's gauge the scalar potential is deter- 
mined by the instantaneous value of the charge distribution, using an equation 
similar to Poisson's expression in electrostatics. The vector potential, however, 
is considerably more difficult to calculate. According (2.19), a time change in p 
implies an instantaneous change in <F This fact denotes the non-physical nature 
of $ since real physical magnitudes can change only after a delay determined 
by the propagation time between the perturbation and the measurement point. 

In this book we will use only the Lorenz condition, but it should be made 
clear that the E and H helds calculated from the potentials with the Coulomb 
or Lorenz gauges must be identical. 

The complete solutions of the inhomogeneous wave equations for the poten- 
tials (2.14) are linear combinations of the particular solutions and of the general 
solutions for the corresponding homogeneous wave equations. The next section 
is devoted to finding these particular solutions, which express the potentials in 
terms of integrals over the source distributions J and p. 

2.2 Solution of the inhomogeneous wave equa- 
tion for potentials 

Let us now calculate the expression of the potentials created by an arbitrary 
bounded source distribution (charges and currents) in an unbounded homoge- 
neous, time-invariant, linear and isotropic medium of conductivity zero that 
we assume to be free space (Fig. 2.1). From (2.14) we see that the scalar po- 
tential <& as well as each of the three components A i7 (i = 1,2,3), of the vector 
potential A satisfy inhomogeneous scalar wave equations with the general form 

□*(r,t) = V 2 *(f,t) - l ^d^ = -9(r,t) (2.21) 

where the operator □ acts on the coordinates r, t of the held point, while the 
sources coordinates are r',t'. 

To facilitate the solution of this equation, we can use, owing to the linearity 
of the problem, the superposition principle and consider a source distribution 



(2.19) 
(2.20) 
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g (r, t) as constructed from a sum of weighted space-time Dirac delta function 
sources, i.e., 

g(f,t)= f f g(f',t')S(f-r')S(t-t')dv'dt' (2.22) 

where V is a volume containing all the sources. Thus, (2.21) can be solved in 
two steps, using Green's method in the time domain, as follows. 

a) The first step is to calculate the response, G(f, r ' , t, t'), generated by the 
space-time Dirac 5— function source, S(f—f')5(t — t'), located at position f' and 
applied at time t' which obeys the inhomogeneous wave equation 

□G(r, r ', t, H) = V 2 G(f, f, t, f ) 1 ^^11 = - § (f- r')5(t U) 

(2.23) 

and satisfies the boundary conditions of the problem. The function G(r, f', t, t') 
is called Green's free-space function, which, because of the homogeneity of the 
space, must be a spherical wave centred at position r' at time t' . This function 
depends on the relative distance, R = \f—r'\, between the point source and 
the observation or field point and on the time difference r = t — t' . Thus 
G(r,r',t,t') = G(R,t) and (2.23) can be written, using spherical coordinates, 
as 

nG(R,r) = ^ d ~^ 1 ~ = S(R)6(r) (2.24) 

where, (??), 

1 d 2 (RG) , nnr . , 

d 2 G d 2 G 

- ~w (2 - 25b) 

b) The second step is to find ^>(f,t) from Green's function. Owing to the 
linearity of the problem and, from (2.22), if the solution of (2.23) is G, then the 
solution of (2.21) is 3 

*= / / g{r',t')G{R,T)dv l dt l . (2.26) 

Jt'=-oo JV 

Because G fulfils the boundary conditions, so too does \I/(r, t). 

To find the Green's function let us consider first a general point R ^ such 
that equation (2.24) simplifies to 

DG{R ' T) - R-^W~ °- (2 ' 27) 

Multiplying this equation by R and defining G' — RG we have the homogeneous 
wave equation 

d 2 G' 1 d 2 G' n 

-dW-^ = ° (2 - 28) 



3 Note that Eq. (2.27) represents the spatial and temporal convolution of g(r,i) and G(r,t) 
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The general solution of the above expression, as can be verified by direct sub- 
stitution, is 

G'(R,t) = f{T-R/c) + h(r + R/c) (2.29) 

where / (t — R/c) and h(r + R/c) are two arbitrary functions of their respective 
arguments and they represent waves propagating along R in the positive and 
negative directions, respectively. Therefore 

G{ R, T) = (±-m + <±±m (2 .3o) 

R R 

The potential that results from substituting Green's function h(r + R/c)/ R in 
(2.26) is termed the advanced potential and is a function of the value of the 
sources at the future observation instant. This advanced potential is clearly 
not consistent with our ideas about causality, according to which the potential 
at (t, r) can depend only on sources at earlier times. Thus, in (2.30) we must 
consider only the retarded /(r — R/c)/R solution as physically meaningful. 

To determine f(r — R/c)/R, we integrate the differential equation (2.23) in 
a very small volume around the singular point R = 0. Thus, taking into account 
that for R — > the function G behaves as f(r)/R, we have 



= - [ S(r-r')S(T)dv' = -6(t) (2.33) 
Jv 

or, since V 2 (1/R) = -4ir5(R) and dv' = AnR 2 dR, 

-J i7rf(r)S(R)dv' + J j R^fidR = -S(t) (2.34) 
As R — > 0, the second integral can be eliminated and therefore 

f(r) = ^ (2.35) 
As the function / depends on r — R/c and /(r) = /(r — R/c)\ R=Q , we have 

f(r-R/c)= 5 ±^ (2.36) 
and the solution of (2.24) is given by 



(2.37) 
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This is Green's time-dependent retarded function, which takes into account the 
time needed for the electromagnetic perturbation to reach the observation point 
from the point source. Substituting this function in (2.26), we have 

W,t) = f J 9 <f J) KT ^R C) dv ' dt '- ( 2 - 38 ) 

and, integrating in t', we finally find that, under the assumption of causality, 
the solution of the inhomogeneous wave equation for potentials is given by 




(2.39) 



where 

[9]=9f,t-*)=g(r',t') (2.40) 

is the value of the source densities evaluated at the retarded times t' = t — R/c, 
which in general are different for each source point, R/c being the delay time 
due to the finite propagation velocity of the electromagnetic perturbations. In 
the following the physical magnitudes evaluated in retarded times are shown in 
brackets. 

By analogy with (2.39) the solutions to the inhomogeneous equations for the 
potentials are 

* Kt > = ii>' < 2 - 41 *> 

*™ - WP <" ib > 

where the bracket symbol [ ] indicates that the enclosed magnitude must be 
evaluated at the retarded time t' = t — R/c. That is 

[ P \ = p(r',t') = p(r',t-R/c) (2.42) 
[J] = J(r',t') = J(f',t-R/c) (2.43) 

are the charge and current densities, respectively, evaluated in the retarded 
times t'. 

Expressions (2.41a) and (2.41b), which are called retarded potentials, in- 
dicate that the potentials created by a distribution at the field point P are 
determined, at a given time t, by the values of the the charge and current den- 
sities at the source points evaluated at previous times t' , which generally differ 
for each source poin. It is easy to check that these potentials, together with the 
continuity equation, verify Lorenz's condition (2.10). 

It should be noted that (2.39) is a particular solution of (2.21), to which a 
complementary solution of the homogeneous wave equation □'^(r, t) = can 
be added in order to arrive at other possible solutions of (2.39). Thus, other 
conditions must be imposed to ensure that the only possible solution of (2.21) 
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is (2.39). These conditions, establishing the uniqueness of (2.39), can be found 
in Appendix ??. 

For sources with time-harmonic dependence 

p(r',t) = Re{p(rV"*)} (2.44) 

J{f',t) = Re{J(f')e jut } (2.45) 

the expressions of the retarded potentials $ and A simplify to 

A(r,t) = ^ Re {/ ^JlVV" (t ~ f Wj ^ Re{A(r)e J "*} (2.46a) 
$(r,t) = — !— Re(/ ^pir'y^-^dv'} = Re {*(f)e^'}(2.46b) 

where 



Mr) = ir- [ ^J(f')e- jkR dv' (2.47a) 

47T Jy, R 

*(f) = — — / 4p(Oe^ fci W (2.47b) 
47re Jyi R 



where k = ui/c = 2ir/\ is the wavenumber in the unbounded medium and A is 
the wavelength in the medium. For harmonic signals the time delay R/c, when 
multiplied by u>, becomes a phase shift given by kR. 

2.3 Electromagnetic fields from a bounded source 
distribution 

The fields created by a bounded source distribution (charges and currents in 
free space) of arbitrary time dependence can be determined by inserting (2.41a) 
and (2.41b) into (2.1) and (2.4). Next, we find the expression for the magnetic 
field first and for the electric field afterwards 4 . 

Magnetic field 

Starting from the equation 

B = V x 1= ^ f V x^-dv' (2.48) 
4tt Jy, R 

and transforming the integrand by the vector analysis formulas (??) and (??) 
of Appendix ??, with \Jj = 1/R and A = J, we can directly find the magnetic 
field equation 

4 An alternative way of obtaining the electromagnetic fields is indicated in Section ??. 
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where [ Jj is the retarded current density at the source point r ' and d Jj dt 

d[J]/dt' — d[J)/dt is its time derivative at the instant t' = t — R/c. 

Expression (2.49) can be written as the sum of the two components B = 
Bbs + B ra d, which are defined below. 

The Biot-Savart term, B os : 



(2.50) 



which is formally analogous to the Biot-Savart expression of magnetostatics, 
although here with the sources evaluated at the retarded times. As this term 
decreases with l/R 2 , its contribution is appreciable only at short distances. 
The radiation term, B ra d'- 



B. 



rad 



4ttc 



x R 



R 2 



-dv' 



(2.51) 



which depends on l/R, and consequently its contribution to the magnetic field 
predominates at long distances from the sources. 

At the static limit, when the sources do not change with time (i.e., for a 
stationary current distribution) equation (2.49) simplifies to the Biot-Savart 
expression of magnetostatics 



= Mo f 

47T Jy, 



J x R 



dv' 



(2.52) 



Electric field 

From (2.4) and (2.41) we see that 



E=-^( V [ 4dv'-^[ §- [ 4dv' (2.53) 
47T£ J v , R 4tt J v , dt R 

Taking into account that d/dt 1 = d/dt and that V* (R) = (d^/dR)VR we 
have 



V R = W V ^ + S VW = [P] 



R 
dR 

v| . [p] 



d [p] dt' 




(2.54a) 
(2.54b) 
(2.54c) 
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which, when substituted in (2.53), and taking into account the continuity equa- 
tion V • J = —dp/dt, gives 

= J- / - lM - 4- [V • J\) dv> (2.55) 

v ' ; 47T£ J v , \ R A c 2 R R 2 c I J J v y 

To get the exact form of the radiation term, which depends on the distance as 
1/R, we need to transform the integrand of this expression by developing V' • [J] 

as 5 

[V • J] = V ■ [J] - ?±M1 (2.56) 
cR 

thus we can rewrite the third term on the right-hand side of (2.55) as 



R 

y,Wc 



V -J 



dv' = - 



R 



V • [J)dv' + 



c 2 K A 



-dv' (2.57) 



The calculation of the first term on the right-hand side can be facilitated by 
calculating just one component, for example the x component 



Rx 



[J] dv' 



c R z 



-I 

-L 



cR 2 



([J]-R)Rx 
cR^ 



(2.58) 



where we have used (??), applied the divergence theorem, and integrated over 
an external surface that encloses the sources in which [jj = 0. Therefore, 
generalizing to three dimensions and inserting the result in (2.55), we get 



4ne, 



( 



+ 



'2([f\- R) R-[J\ (r-r)\ 

v ' v ' ' dr' ■ 



cR 4 



x R) x R 



R 3 



-dv' 



(2.59) 



5 V ' ■ [J] = (V ■ J) t , + ■ V'f = (V ■ J) t , 
Thus 



9[J] v ./ _ rv / fv , , H d[J] 
SV Vt ~ *- V J n' + ~dt r 



(V ■ J) t , = [V ■ J] = V' ■ [J] 



cH 
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which can be expressed as the sum of the three components E = E c + Ei + E r ad, 
which are defined below. 
Coulomb's term, E c , 

This term is similar to the static Coulomb's expression except concerning the 
time delay. 

Induction term, Ei, 



4ire J v , I cR 4 cR? J v ' 

Because of their dependence on 1/R 2 , the contribution to the field of the terms 
(2.60) and (2.61) decrease quickly with distance. 
Radiation term, E ra d, 



R 



Erad wy v , — r 3 — dv -^L 



'21 

dt 



x R) x R 

—dv' 



R 3 

(2.62) 

This term, which depends on 1/R, is the electric field component that predom- 
inates for long distances. Together with (2.51), this component is of interest in 
radiation phenomena (see next subsection) . 

At the static limit, expression (2.59) simplifies to Coulomb's expression of 
electrostatics 

$ = IT— I S dv ' ( 2 - 63 ) 

47T£ J v , R 3 

Alternatively, the electric field can be expressed only in terms of the current 
density, by using the continuity equation. In fact, from (2.56) we have 

[p] = fjV ■ J]dt' = (V • [J] - dt' (2.64) 

Inserting (2.64) into (2.59) and operating in a similar way to (2.58), we obtain 
another alternative expression for the electric field created by a bounded source 
distribution 




(2.65) 
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Fields created by a time-harmonic source distribution 

For time-harmonic dependence of the sources, the field expressions (2.4) and 
(2.1) simplify to 



B = V x A (2.66a) 
E = -V*-jwA (2.66b) 

and equation (2.49) for the magnetic field becomes 

while the different expressions for the electric field, (2.55), (2.59) and (2.65) 
become, respectively, 

1 f p e~ jkR R , . 
E = / i-— dv' + 

47T£ Jy, R 3 

jk f ( P R f(f")\ e~i kR , . 

1 ' • ' dv' (2.68a) 



4tt£ Jy, \ R c j R 

'2{J-R)R 



E 



^s Jy, ^ + 4TCo j v , y cRi cR 2 y + 

ik r (fx R)xR 



(f- R) R 



ATlUJEo Jy, \R 3 R 5 J 

i , (z(j-r)r f\ 

AlTSo Jy, Ci? 4 C i? 2 + 



(jx R) 



+ 7^ / " J e ~ ikRdv ' ( 2 - 68c ) 

A-KEoCjy, R 3 

and the radiation fields (2.51) and (2.62) become 

B = ™±f ^fie-^dv' (2.69a) 

47TC Jy, i? 2 V ' 



r [J x R) x R 

^ / ^ J e~i kR dv' (2.69b) 

4tt Jy, R 3 
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2.3.1 Radiation fields 

Examining the total fields (2.49) and (2.65) generated by a bounded distribution 
of sources with arbitrary time dependence, we find that in general the near-zone 
terms, which depend on 1/R n [n > 1), are negligible compared to the radiation 
terms, (2.51) and (2.62), which depend on 1/R, when the condition 



R » c 



[J] 



d[J}/dt 



(2.70) 



is fulfilled for any of the infinitesimal volume elements into which the source can 
be subdivided. For time-harmonic fields, this condition becomes 

R»X (2.71) 

Hence, the radiation term predominates when distances from the sources 
are great compared to any wave-length involved. The zone where the radiation 
fields predominate can be called by several names: far zone, wave zone and 
Fraunhofer zone. Note that the far zone is farther away from the sources at 
lower time dependence (i.e., at lower frequencies) and there is no far zone at the 
static limit. 

Let us select the reference origin close to or within the source distribution, 
(Fig. 2.1). If the field point is far away from any source point such that r » r' , 
or equivalently r » I, where / is the largest dimension of the source distribu- 
tion, then it is possible to make some general approximations in the expressions 
(2.51) and (2.62) which greatly simplify the calculations. To confirm this, let 
us write R in Fig. 2.1 as 

R=\r-r'\ = (r 2 -2f-f' + r' 2 ) 1 ' 2 (2.72) 

Since the reference origin is close to or within the source distribution, we can 
calculate the radiation fields at distances r » r' by expanding the binomial 
(2.72) as a series in powers of the small parameter r'/r and take only the linear 
terms of the expansion 

1 — 2 — ^ 1 — T ) =r h ... ~ r — r' ■ f = r — r' cosO 

(2.73) 

where 9 is the angle between f and r*. This approximation is equivalent to 
considering that, far away from the sources, r and R become parallel. 

Thus, as r'/r « 1, in the expressions (2.51) and (2.62), we can make the 
approximation 

R~r (2.74) 

in the denominator. This is equivalent to ignoring, in the modulus of the con- 
tribution of each source point to the total field, the difference in the distance 
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travelled by the signal. Thus (2.70) becomes 

[J] 



r >> c 



d[J}/dt 



(2.75) 



and (2.71) becomes 

r » A (2.76) 

In the retarded time, t' = t — R/c, the approximation (2.74) is not valid 
because the sources can be very sensitive to small changes in the delay time 
R/c. Thus, for the delay time, at distances r » r' we need to keep at least 
the two linear terms of the expansion (2.73). Therefore 

. R r f'-f . r' cos 9 .„ „. 

t' = t = t - - + = t' a + 2.77 

c c c c 

where t' = t — r/c. 

Therefore, from (2.77), the retarded time has two components. One, r/c, is 
the time needed for the electromagnetic field to reach the field point from the 
origin of the coordinates. The other, r ' • f/c, represents the time necessary for 
the propagation of the electromagnetic perturbation within the geometric limits 
of the source distribution. This term, given that the largest dimension of the 
source distribution is I, (Fig. 2.1), has a magnitude of 

r" -r/c~l/c«r/c (2.78) 

Hence, using the approximations (2.74) and (2.77) the integrands of the 
radiation fields (2.51) and (2.62) simplify to 

*- - d&Li ' 1 ™*^ (2 - 79b » 

or, for time-harmonic dependence, 

Erad = J (fxrjxr e- jkR dv' (2.80a) 

^ f fx f e-i kR dv' = J ^h [ J x f e-i kR dv' (2.80b) 
47rcr J v , 4-Kcr J v , 

A comparison of Eqs. (2.79a) and (2.79b), shows that the radiation fields 
are perpendicular to each other and to the direction of propagation. They are 
related by 

E = r] H x f (2.81) 
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where the ratio 7] is defined as 

Vo = § =(Mo/£o)^ =120tt n (2.82) 

and is called the intrinsic impedance of free space. 

According to Poynting's theorem the total radiated energy passing through 
the unit area perpendicular to the direction of the vector E rad x H rad is given 

by 



/ Praddt^ f (E rad x H rad )dt (2.83) 

J — oo J — oo 

and the total flow of power passing through the closed surface 5 situated in the 
far-field zone is 



/ / Vrad ■ dsdt = [ [ {E rad x H rad ) ■ dsdt (2.84) 

J-ooJS J—oo J S 

In summary, the assumptions involved in using (2.79) and (2.80) to calculate 
the radiation fields created by a bounded source distribution in the far-field zone 
are: 

a) r » (c\[J]\ I \d\J\jdt\) or, equivalently, r » A for any wavelength of 
the radiation spectrum which allows us to neglect f/r 2 terms. 

b) r » I , where I is the largest dimension of the source distribution which 
allows us to make the approximations (2.74) and (2.77). 



2.3.2 Fields created by an infinitesimal current element 

The simplest case of a bounded source distribution is that of an infinitesimal 
current element i(t), which is assumed to be oriented on the z axis (Fig. 2.3) 
and to have arbitrary time dependence. This current is mathematically defined, 
in terms of the Dirac delta function, as 



— * / \ 7 \7 

J(r, t) = i(t)5(x')6(y')z - — < z' < — (2.85) 

The fields of this current element can be easily calculated by substituting (2.85) 
in (2.49) and (2.65). Thus, we have 
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Fi £Ur6 2 . 3 1 Infinitesimal current element Solo campos de radiacion j jjfalta la r del 
radio vector del punto campo 



Fields created by an infinitesimal current element with arbitrary-time 
dependence: 



H(r,t) 



E{r,t) = 



Az 1 d[i 
47T£„ c 2 r dt 



(f x (f x z)) 



Az 
Att 
Az 
Att 

Az 
4ire 

Az 
4ne n 



1 d[i\ 
or dt 
1 d[i] 
or dt 
1 



+ 



(z x f) 



sm.6 (p 



[i] dt + 



[i] dt + 



Az 1 d[i] 
Atte c 2 r dt 



sint 



(2.86a) 
(3 (z ■ f) f — z) + 
(2cos6» f + smeO ) + 

(2.86b) 



where [i] = i(t — r/c). 



For time-harmonic dependence of the current element, i — Re {le 3 "*}, equa- 
tions (2.86a) and (2.86b) simplify to 
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~7 




Figure 2.4: Radiation field separates from the source and propagates to infinity 
Dibujar el dipolo. Note that there is not radiation in the direction in which the 
current element is pointing. PP 259 del panofsky :-Como puede verse en la figura 
las lineas de campo de raciacion representa una familia de lazos, atravezados 
por las h'neas de campo magnetico, que se propagan hacia el infinito ( i.e. waves 
see chapter tal) 



Fields created by an infinitesimal current element with time- harmonic 
dependence: 



lAz / 1 \ e~ jkr 
4-7T \ jkr J r 

IAz ( , 1 



sin a ip 



(2.87a) 



-jkr 



47T 

lAz 
~2^ 



jkr (krY 



■sin6l + 



■jkVo -7- 



-jkr 



jkr (krf 



■ cos 6 f 



(2.87b) 



These expressions can be also derived directly from the vector potential (2.41b), 
which in this case simplifies to 



A = z^ r ^W^»A: 



(2.88) 



Thus, the magnetic field is given by 6 
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H = — V x 1= — V x (Az) = — (VA x z + A(V x z)) = — VA x z 
Mo Mo Mo Mo 

(2.89) 

where we have applied the vector identity (??) and taken into account that the 
curl of a constant vector is zero. Hence, using spherical coordinates, we get 



H = — VAxz= — — — )rxz 
Hq 4-k or \r J 

Az f 1 d [i] [i] \ „ Az ( \ d [i] 



—— - f x z = —— ; h — sin 6(p 

47r V cr at r z I 47T V cr at r l 



(2.90) 



which of course coincides with (2.86a). 

The electric field (2.86b) can be calculated from (2.90), taking into account 
that from (l.ld), in source-free regions, we have 7 

1 /' 

E(f,t) = — VxH(f,t)dt (2.93) 

£ J-oo 

From the relation between the charge and current, i(t) = dq(t)/dt, we have 

where p = q A z is the dipole moment of a time- varying electric dipole 8 , the 
so-called Hertzian dipole, formed by two point charges with values of +<?(t) 



6 For a given vector field A, the field lines arc defined by the condition that, at any point, 
the line element dl and the field are parallel i.e. dl X A = 0. For the field created by a current 
clement, from Eqs (2.87a). the magnetic field has only (p component and consequently their field 
lines are closed around the Z axis. The radiation electric field has 9 and f components, although 
the radiation electric field has only component which becomes null in the region — > 0. Then in 
this region predominates the E r = E ■ r component and consequently the electric field lines close( 
see Fig 2.4) as would be expected from Maxwell's equations since, outside the sources, there only 
exist curl sources. 

7 Note that once calculcd H = 1/jUqV X A we can obtain E using (l.ld) or (1.67d) and taking 
into account that, in source-free regions, we have 

E=^-JvxHdt = cjv x (V x A) dt (2.91) 

or 

E = — — V x H = —V x (V x A) (2.92) 
jeooj jk 

for arbitrary or harmonic time dependence respectively. Thus wc do not need necessarily to 
calculate $ to obtain the fields. 

The time varying electric dipole is defined as two time varying charges of opposite magnitude 
i<?(t) separated by a constant distance Az much less than the field point r. The dipole moment 
p{t) is given by the magnitude of the charge times the distance A.Z between them and the defined 
direction is toward the positive charge i.e. p(t) = q(t)A.z.. Alternatively it would be possible to 
model the oscillating dipole as two constant point charges of opposite sign separated by oscillat- 
ing distance A.z(t). However, the fields created for such accelerated charges need from the theory 
developed in Chapter ??. 
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and — q(t) and the dot indicates differentiation with respect to time. Thus the 
time- varying current eiement is equivalent to 

= 11 < 2 ' 95 > 

or, for time-harmonic dependence, 

< 2 - 96 > 

Introducing (2.95) into (2.86a) and (2.86b), and (2.96) into (2.87a) and (2.87b), 
we get the field created by an infinitesimal current element (hertzian dipole) in 
terms of its dipole moment as: 

Fields created by a Hertzian dipole with arbitrary-time dependence: 

H = -LfM + M^^ (2.97a) 
Anr \ r c J 

g . i fM + M + M) si „ 9( , + 

4nreo \ r z rc c l J 
Fields created by a Hertzian dipole with time-harmonic dependence: 

h = frG + ^)^- sin ^ (2 - 98a) 
p ( 1 , i k ! 2^ e ~ jkr 



E = -J—l— + J —-k'\ sin 00 + 

47T£ V r T J r 

-^tl J -)— rf (2.98b) 
2ire \r z r J r 

The radiation fields created by an infinitesimal current element can be ex- 
pressed, from (2.86a) to (2.87b), in terms of its current amplitude or of its 
equivalent dipolar moment. 

Radiation fields created by an infinitesimal current element: 

For arbitrary-time dependence 

Az 1 d\i] , 

H rad = 1 -JT Sin 0<P 2 - 99a 

47r cr at 

- slif ^ < 2 - 99b > 
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For time-harmonic dependence 

IAz e~ jkr 

H rad = — — jk sinOtp (2.100a) 

_ TAz e -i kr 

E rad = —jk Vo smee (2.100b) 

47r r 

and from (2.97a) to (2.98b), we have the radiation fields in terms of its equivalent 
Hertzian dipole: 

Radiation fields created by an electric dipole: 

For arbitrary-time dependence 

H rad = J- 00 (2.101a) 
E rad = -r^—H sin 69 (2.101b) 
For time-harmonic dependence 

H rad = f sin 0<p (2.102a) 

47r r 

-» — T)k? e~^ kr 

E rad = -f- sin (96* (2.102b) 

4ne r 

More details about this elemental radiators and how they can be physical 
approximated are given in subsubsection ?? of chapter ??. 

2.3.3 Far-zone approximations for the potentials 

The general expressions (2.49) and (2.59) for the fields due to an arbitrary 
source distribution of finite size are of theoretical and sometimes of practical 
interest. However, except for the case of the infinitesimal current element, it 
is much easier to calculate the fields created by a given source distribution via 
the potentials, as indicated in Fig. 2.2. This can be seen simply by comparing 
the complexity of the expressions for these fields, (2.49) and (2.59), with those 
for the potentials (2.41a) and (2.41b). Because of the vector product in the 
integrand of (2.51) and (2.62), this argument continues being true even when 
we are interested only in the radiation fields. In the far zone, we can make the 
approximations (2.74) and (2.77) for the potentials. Hence, the integrands of 
the retarded potentials (2.41) simplify to 

$(r,t) = -L / P^ldv'c^-^—f p(f', t ' + —)dv' 
4tt£ J v , R 47T£ r J v r yu c 

(2.103a) 
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The magnetic field can now be calculated from (2.1), using (??), as 

4tt./ v , r 4tt./ v , c r 

(2.104) 

where, if we are interested only in the radiation field, the second term can be 
ignored since it depends on 1/r 2 , and therefore 



i /■ 

4tt y v , r 



Mo 47r ./v 

Furthermore, from (??), we have Vx J(*) = V^xdJ/d^ with * = t' +f'-f/c. 
Thus, it follows that 

C COT 

f aj(r',id + ^) 
= — x — ^ 2.106 

cot v y 

and therefore 

H = rx^— (2.107) 

fl C ot 

which, as would be expected, leads to (2.79a). If the time variations of the 
sources are harmonic the expressions (2.103a) , (2.103b) and (2.107) become 

* = —!—e-3 kr [ p{r')e^ r ' dv' (2.108a) 

47T£ r Jy, 

A = £L e -Jkr J j{f') e ok-f' dv , (2.108b) 
H = -— f x A (2.108c) 

The radiation electric field can be calculated from (2.107) or (2.108c) simply 
using (2.81). 



2.4 Multipole expansion for potentials 

In many cases, such as the study of most antennas, in order to calculate the 
radiation fields, we cannot make any approximation concerning the potentials 
other than those assumed above. For example, we need to carry out the integra- 
tion in (2.103b) or (2.108b) in order to calculate the vector potential. However, 
if we assume that the charge distribution does not change appreciably over time 
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l/c, we may expand the integrands of (2.103) in a Taylor series about t' in terms 
of the parameter f ■ r/c. For example, for the vector potential, we have 



r ■ r 



+ ... 



(2.109) 



where we have omitted higher-order terms in f'-f/c. Thus after inserting (2.109) 
in (2.103b), we can write A as the power-series expansion 



A 1 + A 2 + ... 

4irr 



/i0 / J(f',t' )dv'+ "° 
Jv 



Ancr 



dJ(f',t') 



at' 



f 1 ■ fdv 1 + ... 

(2.110) 



Therefore the first two terms of the expansion (2.110), A\ and A 2 , are given by 

j v J(f>,Qdv' (2.111a) 

Mo f df(r',t> ) 
4-ircr J v , dt' 



Ai = 

M = 



f' ■ fdv 1 



Mo 



(2.111b) 



t'=t' 



If the time dependence of the sources is sinusoidal the condition that the source 
distribution does not change appreciably over time I /c is equivalent to assuming 
that l/c « T (where T is the period of the signal) or equivalently l/X « 1, 
i.e., that the dimension of wavelength is much greater than that of the source 
distribution 

A »Z 



In this case, we can perform the series expansion 



Jk-r' _ jkf-r 1 



1 



1+jkr ■ r' - -k 2 (f ■ f') 2 + ... 



which, after substituting in (2.108b), leads to 

A = A i + A 2 + 

where 



A x 
A 2 



Mo e 

47T 



-jkr {• 

— \ JVW 

r Jv 1 



= jk 



Mo e " 



-jkr 



An r 



[ J{r')r ■ f'dv' 
Jv 



(2.112) 
(2.113) 

(2.114) 

(2.115a) 
(2.115b) 
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which are the Fourier transforms of (2.111a) and (2.111b), respectively. 
Of course, there are analogous expressions for the terms of <& 

<& = $1 + $2 + •■■ 



Note that, since the contribution of each point source to the integral in 
(2.117a) is evaluated at the same time t' , this integral represents the total 
charge of the source distribution. Thus, if the net charge of the distribution is 
zero, we have $i = 0. If the net charge is not zero, the constant, the electrostatic 
potential $i created by that charge depends on r~ 2 and consequently it does 
not contribute to the radiation. 

The expansion (2.110) allows us to decompose the electromagnetic field 
created by a time-varying source distribution of finite dimension in terms of 
elementary time- varying source distributions, called electric and magnetic mul- 
tipoles, located at the origin. This is similar to the well-known multipolar ex- 
pansion of the electrostatics (or magnetostatics) to decompose the field created 
by a stationary source distribution of charge (or current) in terms of electric (or 
magnetic) multipoles. However, now the original distribution is time-varying 
and produces both electric and magnetic fields. Thus, as result of the expan- 
sion, we will obtain both, electric and magnetic multipoles. To verify this, we 
next analyze the first two terms, (2.111a) and (2.111b), of (2.110). 

2.4.1 Electric dipolar radiation 

The evaluation of the term (2.111a) of the power-series expansion of A can be 
facilitated by calculating just one component of J v , J(f' , t' Q )dv', for example the 
x component 




(2.116) 



where 




(2.117b) 



(2.117a) 





(2.118) 



since 




(2.119) 
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as can be seen by applying the divergence theorem and by integrating over 

an external surface, where J(f',t' ) = 0, that encloses the sources. Therefore, 
generalizing to three dimensions we have 

/ J(f',t' )dv' = - [ rV • J(f' ,t' )dv' (2.120) 
Jv Jv 

and using the equation of continuity 

V'-J(r', i ' ) = -^J^) (2 . 121) 



we get 



iv Jv 
which, when substituted in (2.111a), gives 



J v J(f',t' )dv' = j^r' ^f^ dv' (2.122) 



M = ^( f' ^'^ dv' = -f* £ [ f'p(f',t' )dv' (2.123) 
Attt J v , dt Airrdt J v , FK ' °' V ' 

The integral J v , r' p{f' ,t' Q )dv' is by definition the electric dipole moment, [p\, 
evaluated at the retarded time t' , of the time- varying source distribution, i.e., 

[p\ = [ f"p(f',t' )dv' = [ f'p(r',t--)dv' (2.124) 
Jv' Jv c 



Thus we have 



A p L m = ^P\ (2125) 

4vrr dt 4?rr v ; 



The magnetic radiation field, from (2.107), is given by 



tlrad = -—. = — j <P (Z.tZb) 

where we have assumed the direction of p parallel to the polar z axis. This 
expression, as might be expected, coincides with the radiation term, (2.101a), 
of (2.97a). From (2.126), the electric radiation field, given by (2.101b), can 
be obtained using (2.81). Of course the corresponding expressions for time- 
harmonic fields are given by (2.102a) and (2.102b). Therefore, in a preliminary 
approximation, the original source distribution can be replaced by an electric 
dipole located at the origin of coordinates. 



2.4. MULTIPOLE EXPANSION FOR POTENTIALS 



55 



2.4.2 Magnetic dipolar radiation 

The analysis of the term (2.111b), can be facilitated by expressing the integrand 
as follows 

J(f',t' )(r.f') = \(j(f',if )(f'-f)-f'(j(f',t' ).r)) 
+ \ (J(r\ t'oW ■ f) + f (J(r ', f ) • r)) 

= \r x ( J(r ', t' ) x r ') + \ ( J(r ', t' )(f ■ r>) + r ' ( J(f ', ^ • 

(2. 

Then, substituting in we get 

^2 = A 2 m+A 2g (2.128) 

where 



4 ^° 
and 



A ' 29 = 8^1 / v , ( J V^o)(^ • O + r' (j(f'X) • r)) dv' (2.130) 



The integral (2.129) can be written as 



i2m 



47rcr 9i 
where 

\m 



x f (2.131) 



= / ^f^W (2.132) 

is by dehnition the magnetic dipolar moment about O, evaluated at the retarded 
time t' 0> of the source distribution. Thus, under the assumption that m = mz, 
the magnetic radiation field given by (2.107) is 

H rad =-L-f x(f x[m]) = -!-Msin^ (2.133) 
From (2.81) the electric radiation field is given by 

E rad = --^- [ ^sme(p (2.134) 
Attt c 

For time-harmonic dependence, we have 

H rad = d^i e -i^ (2.135) 
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and 

-. fc 2 msin# _,- fer „ , n -,nn\ 
E rad = — -. Vo e J <P ( 2 -!36) 



where 

m 



These expressions are similar to (2.101a)-(2.102b), which were obtained for 
the electric field of the radiation of the electric dipole. In fact, as we will see 
in the next section, there exists a duality in the analysis of the electric and 
magnetic dipoles. 

In the particular case of a current loop of radius a , Fig. ??, for which the 
current i does not change appreciably over time a/c (or equivalently a << A for 
any frequency involved), (2.132), becomes 

m = i/^=i5 (2.138) 



r 

where T is the countour of loop and S is the vector area of the surface subtended 
by the contour T. In this expression, Jdv has been changed to idl. The surface 
vector S is directed normal to the loop according to the right-hand rule for 
the direction of the current in the loop. Thus, for the circular current loop, the 
radiation fields (2.133)-(2.136), can be written, for arbitrary time dependence, 
as 

H rad = /-^J sin 00 (2.139a) 

E rad = -^Msin^ (2.139b) 

where i is evaluated at t' . These equations, for time- harmonic dependence 
become 

H rad = - fc2 i gsin V ^0 (2.140a) 

471T 

_ k 2 ISsin8 ... x 
Erad = —^- r Vo^ 3 <P (2.140b) 



It should be mentioned that the magnetic moment is important only when 
there exists no radiation of the electric moment of the system. Otherwise the 
one due to the magnetic moment may be ignored. Effectively, comparing Eqs. 
(2.101b) and (2.134), and using E p and E m to indicate the amplitudes of the 
electric radiation fields from an electric and a magnetic dipole, respectively, we 
have 

E Vrn , cp 
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Figure 2.5: Poncr m = iS , poner i en vez de I y el contour I\ A circular loop of 

current in the x-y plane y dibujar campos como en el clcmcnto de corriente. Hacer el dibujo igual 
que el del electrico 



or for time-harmonic variation with both dipoles oscillating at the same fre- 
quency, 



E m rad m 

Since from (2.137) we have 



(2.142) 



m = / T ' * — do' = \ f Po f x udv' (2.143a) 

JV' ^ 2 Jyl 

po = I p^dv 1 (2.143b) 
Jv 

and consequently 

mo ~ up (2.144) 

where u is the velocity of motion of the charges. Thus from (2.142) we have, for 
u « c, 

E Prad » E mrad (2.145) 

i.e., the magnetic dipolar radiation may be ignored in comparison with the 
electric dipolar radiation. 

2.4.3 Electric quadrupole radiation 

The second term, A2 q , of A2 in (2.130), is associated with the electric quadru- 
pole radiation, but to see this we must transform it further. To this end let us 
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consider the x component of the first summand of the integral J v , ^ J(r ', t' )f ■ f'^j dv' 



i.e. 



J f-f' (j(r ', *{,) • x') dv' = J f-f ( J(f t' ) • VV) d«' 

= y V ■ (V (f • f ') J(f, <,)) dv' - J x'V • ((r • r ') J(r O) dv' 

(2.146) 

where the first integral is null, as can be seen using the divergence theorem to 
convert the volume integral in a surface integral with the surface of integration 
outside of the source distribution. Thus 

J f-f (j(f t' Q ) ■ f) dv' = -J x'V • ((f • f ') J(f, t£)) dv' 

= - [ x'V(f-f) ■ J(f',t' )dv' 
Jv 

- f x' (f-f')V • J(f',t' )dv' 
Jv 

(2.147) 



but 



therefore 



V'(f-r') = f (2.148a) 
V'-J(rV ) = -^J^ (2T48b) 



/ X(f',^)f-fW 

= - y x' (j{f',i! )-f)dv' + J x' {f ■ f')^^-dv' (2.149) 



iv 

Generalizing to three dimensions 



J J(f, t' ) (f ■ f) dv' = - J f (j(f ', 4) • f) dv'+ ^ f' (f • r ') 
and therefore, substituting in (2.130), we have 



at 

(2.150) 



;)2 

/v 

The magnetic radiation field, given by (2.107), is 
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The above expression can be written in a more useful form by adding the 
term fr' 2 p(f',t' ) to the integrand 

H2 9rad = -^- r r x (3? (f ■ r>) rr' 2 ) p(r>, t> )dv> (2.153) 

Note that, since f x rr' 2 = 0, the added term do no affect to the value of the 
integral. The advantage of including this term is that, now, the integrand can 
be written as the product of a second rank tensor Q, called electric quadrupole- 
moment tensor of the source distribution, and the vector f 

[ (3r" (f • f ') - fr n ) p(f', t' )dv' = [Q]f (2.154) 
Jv 

The elements of [Q] are 

[Q a p]= I (3x' a x' f3 -r' 2 6 a p)p(r',t' )dv' (2.155) 
Jv 

and [Q]f is a vector with components 

£[Qa/j]r/s (2.156) 

a 

Therefore the radiation magnetic field from a varying electric quadrupole is 
given by 

or, for, time-harmonic dependence, 

H2 Qrad = ^-/^- kr) r x Qf (2.158) 

The radiation electric field can be calculated as usual by (2.81). ft can be 
shown that quadrupole radiation fields are of the same order as the magnetic 
dipole moment and thus much less than that corresponding to the Hertzian 
dipole (Ejercicio).. 

Of course, if we continued analyzing other terms in the expansion tal, we 
would find other multipole moments, such as magnetic quadrupole radiation, 
electric octupole radiation, etc. However, for this, other more complex mathe- 
matical methods provide the results more systematically. 



2.5 Maxwell's symmetric equations 

It can be observed from (l.la)-(l.ld) that Maxwell's equations present a certain 
symmetry that, except in free space and with no source terms, is not complete 
because of the absence of magnetic charges and currents. Indeed, despite many 
experimental attempts, no free magnetic charges or monopoles have been found 
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in nature nor, therefore, would magnetic currents be created 9 . Nevertheless, 
from a purely theoretical standpoint, nothing prevents us from assuming the 
existence of magnetic monopoles; therefore, to complete Maxwell's equations 
we must add the necessary magnetic source terms in order to achieve com- 
plete symmetry between electric and magnetic quantities. To this end, we can 
reformulate Faraday's law (1.1c) and Gauss' law for magnetic fields (1.1b) by in- 
troducing, on their right-hand side, hypothetical magnetic current densities J m 
(Vm~ 2 ) and magnetic charge densities p m (Wb/m 3 ), respectively, as additional 
source terms. With these new quantities included, we can rewrite Maxwell's 
equations for the case that both, electric as well as magnetic sources, exist in 
free space, in the following completely symmetric manner: 
Differential form of Maxwell's symmetric equations 

V-f3 = p (2.159a) 
VB = Pm (2.159b) 

Vx£ = -J m -Po-^- (2.159c) 
dE 

VxH = J + e — (2.159d) 
at 

Integral form of Maxwell's symmetric equations 

D-ds = Q T (2.160a) 

s 

B-ds = Q m (2.160b) 

s 

j> E-dl = - J J m -ds-^ J B-ds (2.160c) 

j> H-dl = j J-ds+^ j D-ds (2.160d) 

It should be emphasized that the symmetrization of Maxwell's equations 
is a powerful mathematical tool which greatly facilitates the solution of many 
practical problems such as the radiation and scattering from aperture antennas 
or permeable bodies. 

Taking the divergence of (2.159c) and using (2.159b) 

V ■ V x E = -V ■ J m - 9V gt B = (2.161) 



^It should be emphasized that, although there is no experimental evidence for the existence of 
magnetic charges, such existence does not violate any known principle of physics. In fact, from a 
purely theoretical viewpoint, Dirac showed [P. A.M. Dirac, Proc Roy. Soc.Lond. A133, 60 (1931)] 
that the existence of magnetic monopoles with magnetic charge g would explain the quantization 
of the electric charge e. We refer to the magnetically charged particles as magnetic monopoles or 
simply monopoles. 
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we get the equation of continuity 

V-J m = -^ (2.162) 

which expresses the conservation of magnetic monopoles and has the same form 
as that for the electric charges (1.3). 

In linear media, we can apply the superposition principle and split each one 
of the field quantities, E, D, H and B, into the sum of two components 



D = D e + D m = e (E e + £ m ) = e Q E (2.163a) 
B = B e + B m = n [S e + # ro ) = p H (2.163b) 

where the quantities with the e subscript depend only on the "true" electric 
sources p and J while the quantities with the m subscript depend only on the 
"hypothetical" magnetic sources p m and J m . In this way, we divide Maxwell's 
equations into two groups corresponding to the field components associated with 
the electrical and magnetic sources, respectively; that is 

V-D e = p (2.164a) 

V-B e = (2.164b) 
Off 

VxE e = -Vo-gf (2.164c) 
BE 

VxJJ e = J + eo-gf (2.164d) 



V-D m = (2.165a) 
V-B m = Pm (2.165b) 

f)ff 

Vx£ m = -J ro - Mo -^ (2.165c) 
BE 

VxH m = e -^ (2.165d) 

Note that the sum of each expression (2.164), added to its equivalent (2.165), 
gives (2.159) and that the set (2.164) coincides with the conventional Maxwell's 
equations (2.159), and that Eqs. (2.164) are formally identical to Eqs. (1.1a)- 
(l.ld) and therefore can be solved as in the previous sections by means of the 
scalar and vector potentials $ and A. Thus, from (2.4) and (2.1), we have 

B e = V x A (2.166) 
BA 

E e = -V$--^ (2.167) 
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where 



V-A + n e — =0 



and where A and $ fulfil the wave equations (2.14a) and (2.14b) 

d 2 A 



V 2 X 



Po £ o 



V 2 $ - Mo^o 



dt 2 
W 



-H J 

P_ 

£o 



the solutions to which are the retarded potentials (2.41a) and (2.41b) 

i r [p\ 



$ = 

A = 



47T£ Jy R 

Mo f [J] 



47T 



R 



dv' 
dv' 



(2.168) 

(2.169) 
(2.170) 

(2.171) 
(2.172) 



The fields created by the magnetic sources p m and J m can be deduced by ob- 
serving that equations (2.164) are transformed into (2.165) and vice versa with 
the simultaneous replacement of the following quantities, called duals 



(2.173) 



The fields E e and H e associated with the electric sources can be calculated 
from the magnetic vector potential A and the electric scalar potential $ by 
means of (2.171) and (2.172). To calculate the fields H m and E m we can use the 
same formalism defining two new potentials, termed "electric vector potential" 
F and "magnetic scalar potential" ip, such that 



Ee 


dual of 


H m 


He 


dual of 


Em 


£o 


dual of 


Mo 


Mo 


dual of 


£o 


P 


dual of 


Prn 


J 


dual of 


Jm 



A 



dual of F 
dual of ip 



Hence 



F = 



^Mo Jv> 

L 



[Pr, 



47T W 

£0_ f [Jm 

4ttL R 



R 
] -dv' 



dv' 



(2.174) 

(2.175) 
(2.176) 



which are the dual expressions of (2.171) and (2.172). 
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By substituting the magnitudes in the first column of (2.f73 and 2.174) for 
their duals in the equations from (2.166) to (2.172) we get 

D m = £A = -VxF (2.177a) 
dF 

H m = -VV-^- (2.177b) 

V-F + Mo £o^ = (2.177c) 

in which if) and F satisfy wave equations that are analogous to (2.169) and 
(2.170): 

d 2 F 

V 2 F-Mo^o^ = -ej m (2.178) 

Thus, by the superposition principle, if both current densities J and J m exist 
simultaneously in a region of free space, the total field E produced at any point 
is the sum of E e and E m given by (2.167) and (2.177a). Hence 

-> dA 1 

E = E e +E m = -V$- — VxF = 

at e 

BA 1 

V • Adt — V x F (2.180) 

dt e 



where Lorenz gauge Eq. (2.10) has been used to express E in terms of A and 
F. 

The total field H is determined analogously from (2.166) and (2.177b) 

- dF 1 - 

H = H e +H m = -V^- — + —X7xA (2.181) 

dt fi 

In practice, it is not necessary to use the latter expression, because once E 
has been calculated using (2.180), by substituting the result in (2.159c), with 
J m = we obtain H. 

2.5.1 Boundary conditions 

It is easy to show, ejercicio, that the boundary conditions corresponding to Maxwell's 
symmetric equations are a logical extension of (1.35); that is, 

n 



(pi-Di) = P s (2.182a) 

n-(Bi-B 2 ) = Psm (2.182b) 

h x (Et - E 2 ^j = -J sm (2.182c) 

(iTi-# 2 ) = J s (2.182d) 



n x 
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in which h is the normal unit vector that goes from region 2 to region 1. Equa- 
tions (2.182b) and (2.182c) show the additional effects of the imaginary sur- 
face magnetic charges and currents, p sm and J sm , at the interface. According 
to (2.182c) and (2.182d), the tangential components of the fields on a real or 
imaginary surface S can be written in terms of surface distributions of electric 
currents 



= J s (2.183) 
- J sm (2.184) 



fix H 

and magnetic ones 

- fix E 

s 

2.5.2 Harmonic variations 

For harmonic variations, the symmetric equations (2.159) simplify to 

V D = p (2.185a) 

VB = p m (2.185b) 

VxE = -f m -jfj, ujH (2.185c) 

VxH = J+jeouE (2.185d) 

and the wave equations for the magnetic scalar potential tp, the electric vector 
potential F, and the Lorenz relations are 

V^ + wVo^ = -— (2.186a) 

Mo 

V 2 F +u) 2 fi Q e F = -£ Jm (2.186b) 
?'V • F 

ip = J - (2.186c) 

weoMo 

with the solutions to (2.186a) and (2.186b) being 

F = £i / dv > (2.188) 

4tt J v , R 

The total field E produced at any point is the sum of E e and E m , and is 
given by 



E = E e + E m = —j — V ( V • A ) - juA V x F (2.189) 

UJ V / £n 



,2 

uj \ :,) 
while for the total field H we have 



H = — j — V f V • — jcoF + —V x A. (2.190) 

uj V / fj, 

where A is given by (2.47a). 
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2.5.3 Fields created by an infinitesimal magnetic current 
element 

From (2.86a) and (2.86b), using the dual equations (2.173), we deduce that the 
fields generated by an infinitesimal magnetic current element, 

_ A 7 A 7 

J m {f, t) = i m (t)S(x')S(y')z - — <z '<— (2.191) 

are given by, (Fig. 2.6), 
-. Az ( 1 d [i m ] [i m ] \ . „ 

E = -TT y-^ + hr )smO<P (2.192a) 



47T \cr rfi r 2 

a - ^(^/j iJ<i(+ ^) (2 » 8,,f+si " w)+ 3l:i^ 8i " w 

(2.192b) 

or, for time-harmonic variation 

E = _^I^-fc(i + £lf!lsin0 (2.193a) 
47r \ jfcr J r 

_ I m Az / 1 1 \ e - - 770 '" - I m Az / 1 1 \ e -i fcr 

47T ^ W£ ° y ^ jkr k 2 r 2 J r Sm 27r ■-' weo ^ ^2 r 2^ j^J r 

(2.193b) 

Comparing the radiation terms of these equations to (2.139a)-(2.140b), we find 
that 

di 

i m Az = n S- (2.194) 
or for time-harmonic dependence. 

I m Az = jujfj, Q IS (2.195) 



2.6 Theorem of uniqueness 

Whenever we have to resolve a differential equation, it is desirable to know 
the conditions that must be fulfilled in order to state that a unique solution is 
possible. In our context, this means to seek the conditions for which we can state 
that there exists a single electromagnetic field that satisfies, simultaneously, 
Maxwell's equations and the given boundary conditions. 

Next, we establish these conditions for non- harmonic and time-harmonic 
electromagnetic fields. 
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X 




Figure 2.6: solo los campos de radiacion se representan" ' 
2.6.1 Non- harmonic electromagnetic field 

A non-harmonic electromagnetic field that varies in a linear region V bounded 
by a surface S is uniquely determined from an initial time, t = to, if the following 
are known: 

i) The values of the sources at each point and at each time for every t > t 
within the region. 

ii) The values of the electromagnetic field (E and H) at each point of V at 
the initial time t — to- 

iii) The tangential components of the electric field E or of the magnetic field 
H on the entire the surface S for all t > to, or, alternatively, the tangential 
components of the electric field E in any part of S and of the magnetic field 
H in the remaining part of S, for all t > to- 

Proof 

This theorem can be proven by a reduction to absurdity- that is, by showing 
that to assume the opposite of what is postulated would lead to a contradiction. 
Let us assume that having defined the three above conditions within a volume 
V, there exist two different electromagnetic fields, (E\ and Hi) and [E 2 and 
H 2 ), respectively, which are solutions to the problem. Given the linearity of 
Maxwell's equations, any linear combination of these two solutions must in itself 
be a solution. In particular, the difference between the two aforementioned 
fields, i.e. the field defined by (& = E 1 - E 2 and H' — H 1 — H 2 ), must also 
be a solution to the problem. Given that, from the hypothesis, the sources are 
the same for the fields {E\ and H\) and {E 2 and J7 2 ), the field {E',H') is 
source-free in V. Thus, if we apply the Poynting theorem (1.39) to (E',H'), we 
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get 

° = it J 'D' + B' ■ H')dv + J aE' 2 dv + j>{E' x H') ■ ds (2.196) 

It is straightforward to show that if the tangential components of the electric 
field E and/or of the magnetic field H are uniquely determined on surface S, 
the final term in (2.196) is null. By integrating this expression with respect to 
the time from t to t and, taking into account that the initial values for t = t 
are defined for all V, we find that 

= J ^(E' ■ D' + B' ■ H')dv + (^j aE' 2 dv^dt (2.197) 

As both of the terms on the second member in (2.197) are positive, this 
equality can be fulfilled only when both E' and H' are null (i.e. when E\ = E 2 
and Hi = Hi), which is what we set out to prove. 

2.6.2 Time-harmonic fields 

In the case of harmonic variations, the uniqueness theorem states that a field 
in a lossy (a ^ 0) 10 region is uniquely determined by the sources within the 
region together with the tangential components of the electric field E or of 
the magnetic field H on S, or, alternatively, the tangential components of the 
electric field E in any part of S and of the magnetic field H in the remaining 
part of S. 

Proof By a reasoning similar to that used for the above case, but using the 
expression (1.111), we get 

By making the real and the imaginary parts equal to zero, we see that these 
two equalities imply that H' and E' are both equal to zero only if a ^ 0. This 
is why we started from the premise that the medium occupying the volume has 
a conductivity that may be arbitrarily small but which is non-zero at all points. 
The field in a lossless region can be considered the limit to the lossy case when 
such losses tend to zero. 



10 The reason why we need the extra condition of the space to be lossy for time-harmonic 
signals is that, by definition, a pure harmonic signal has an infinite duration. 
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Chapter 3 

??Electromagnetic waves 



In chapter 2 the fields created by a bounded time-variyng source distribution 
were calculated and in particular we found that the radiation field propagates 
energy far away from the sources. Of all the possible solutions for the wave 
equation, we will examine primarily the properties of their plane-wave solutions, 
i.e. , waves for which the wave-front are planes 1 . Plane waves constitute a good 
approximation to actual waves in many situations because at sufficiently large 
distances from the sources, in a sufficiently small region, any wave front can be 
treated as a plane wave. For example, a great deal of optics is founded on the 
plane-wave approximation and, similarly, in radiocommunications the radiated 
field at sufficient distance from the antenna can be considered to be a plane 

wave. Moreover, it is possible to demostrate that, in general, an electromagnetic field can puede 
descomponerse como snma lineal de ondas planas { see Appendix ?? ) In this Chapter we consider 
this kind of waves in a linear homogeneous isotropic medium librc de fuentcs. Then incidencia 
normal y oblicua. Ondas esfericas , desarrollo en ondas planas? 

Harmonic. .Electromagnetic waves are not limited in wavelength and in fact cover the spec- 
trum from gamma rays (wavelengths of an n 10-12 cm???????) through X-rays, visible light, 
microwaves, and radio waves, to long waves (hundreds of kilometers long). 

3.1 Wave equation 

For time- varying electromagnetic fields it is possible to combine Maxwell's equa- 
tions to eliminate one of the fields, H or E, to obtain two uncoupled second-order 
differential equations, one in E and the other in H, known as wave equations. 
To formulate these wave equations, let us consider a non-magnetic (p = /x ), 
homogeneous, linear and isotropic region where, in general, source terms J and 
p may exist. Taking the curl of (1.1c) and using the vector relation (??) we 



1 Wave- front is defined as a surface that, at any time t, is orthogonal to the propagation 
vector h at all the points on the surface. 
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have 



VxVx£ = V(V • E) - V 2 E = <A " b 



dt 



9,7 ? 9D, 



- 3 -> 3D 
V 2 E = V(V-E) + ^-(J c + J+—) 

1 dE df d 2 E 



(3.1) 



where J and J c are the source and induced conduction density of the currents, 
respectively. Thus, rearranging terms, we get 

VE -Wto-w-ur = + (3-2) 

which is known as the inhomogeneous vector- wave equation for the electric field. 

A similar equation can be written for the magnetic field H by taking the 
curl of (l.ld), 

Off 0^ ff 
V 2 i? - ^o— - Mo£-^r = -V x J (3.3) 

For a lossless media (3.2) and (3.3) reduce to 

0^ E 1 3 J 

V 2 ^-m £^ = -VxJ (3.4b) 

These Eqs are analogous to the inhomogeneous wave equation for the vector 
potential (2.14a), and consequently their solutions take the form of the retarded 
vector potential given by Eq. (2.41b), i.e. 



E{r,t) 
M(r,t) 



i r + 



— / 

47T£ Jy, 
-l 

^ Jv> 



dv' 



J 



R 



(3.5a) 
(3.5b) 



from which, by means of straightforward operations, we can obtain the expres- 
sions (2.49) and (2.55) for the fields created by a bounded distribution of finite 
densities of charges and currents with arbitrary space and time dependence. 

In source-free regions ( J = 0; p = 0, except the charge and current densi- 
ties induced by the presence of the fields, which are expressed in terms of the 
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constitutive parameters) the equations (3.2) and (3.3) simplify to 

d 2 E dE 
V2 ^-^^r-Mo^ = (3.6a) 

&^ H dH 
V 2 H - u, e— - u, a— = (3.6b) 

which are the homogeneous wave equations that determine the propagation of 
the fields E and if in a sourceless homogeneous, linear and isotropic medium. 
The solutions to these wave equations must be compatible with Maxwell's equa- 
tions and the coefficients of the solutions must be derived from the boundary 
conditions. 

Uniform plane waves are defined as waves with a field amplitude that, at 
any instant, is the same at all points of the wave-front plane. Thus, the field 
amplitude depends only on the distance £ from the origin to the plane (fig. 6.1). 
Therefore, if n — £/£ is the unit vector that is normal to the plane, the del 
operator V becomes V = d/d£ h and Maxwell's equations simplify to 

h~ = (3.7a) 

n-°— = (3.7b) 

„ dE dB . „ , 

„ dH - dD , nrrl . 

nx- = *E+- (3.7d) 

and the wave equations become 

(3.8a) 
(3.8b) 

These equations, which describe the propagation of plane waves in a homoge- 
neous conducting medium, are called the "telegrapher's equations". For nondis- 
sipative media, for example the free space, these equations simplify to 



d 2 E 


d 2 E 


dE 




M(,£ dt 2 




d 2 H 


d 2 H 


dH 


&e 


- ^ dt 2 





d 2 E 


d 2 E 


d 2 E 


1 d 2 E 


of 


M ° £ ° dt 2 


df 


c 2 dt 2 


d 2 H 


d 2 H 


d 2 H 


1 d 2 H 


df 


M ° £ ° dt 2 


df 


c 2 dt 2 



= (3.9a) 
= (3.9b) 
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Figure 3.1: The wave tal, in a lossless medium, propagates at velocity ted to 
the right without changing shape 




Figure 3.2: poner plane wave front 
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3.2 Harmonic waves 

For time-harmonic fields, when the medium presents a conductivity a and, at 
the operating frequency, a complex dielectric constant, e c = e' — je" ,(1.71), the 
wave equation (3.6a) can be written as a time-independent wave equation 

V 2 E - ju/J, o-E + fi uj 2 £ c E = V 2 E - jujfi a e E + fi oj 2 e'E 

= V 2 E +uj 2 u. s / (I -jtm\5 d )E 
= (V 2 + w 2 Mo e ec )E = (3.10) 

where a e = a + ue" ', tan^d = <J e /u>£' ', and e ec = e'(l — j tan<5<j), are the effective 
conductivity, the loss tangent and the effective complex permittivity defined 
in (1.78), (1.81), and (1.83) respectively. 

According to Subsection (??), depending on the characteristics of the medium, 
the values of the term jtan^ in Eq. (3.10) may range from << 1 (zero for a 
perfect dielectric or lossless medium) to >> 1 (infinite for a perfect conductor). 
In a highly conductive medium tan<5d >> 1 and 1 — jtanSd ~ — j tanJ^, and 
thus Eq. (3.10) becomes the so-called time-independent diffusion equation for 
the electric field E 

V 2 E - jujfi aE = (3.11) 

which is of the same type as the one that determines the propagation of heat by 
conduction or by diffusion. As commented in Subsection (??), for most metals 
the relaxation time r is 10~ 14 s, which is a low value compared with the period 
for all frequencies lower than the optical ones. Thus, since t&nSd = (tcj)^ 1 , the 
diffusion equation is adequate for metals at all these frequencies. 
Equation (3.10) can be written more concisely as 

V 2 E- 7 2 E=0 (3.12) 

where 7 is in general a complex quantity called the complex propagation con- 
stant, which, from (3.10) and (3.12), is given by 

-7 = u fj, (s c -j-) 

id 

= w 2 Mo ( £ '-j(e" + -)) 

= uj 2 ^ s' (1 - j tan S d ) = k 2 (1 - j tan 5 d ) = cj 2 fi e ec 

(3.13) 

where 

k = uj^ ! (3.14) 

is the wavenumber corresponding to an unbounded lossless medium with a real 
dielectric constant e' . 

Analogously, for the magnetic field, we have 

V 2 H- 7 2 H = (3.15) 
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3.2.1 Uniform plane harmonic waves 

For uniform plane waves, we have V 2 = <9 2 /<9 2 £ and Eqs. (3.12) and (3.15) 
simplify to 



8 2 E 

W 
d 2 H 



7 2 E = 



7 2 H 



The complex propagation constant 7 is usually written as 2 
7 = jk(l - jtm\5 d ) 1/2 =a+jP 



(3.16a) 
(3.16b) 

(3.17) 



where the imaginary part, [3, is termed the phase constant, whereas the real 
part, a, is called the attenuation constant of the wave. Thus, from, (3.13) and 
(3.17), we can easily calculate the explicit expressions for [3 and a 



(3 = oj 



H e 



i\ 2 



(l+tan 2 5 d ) 1/2 + l 



1/2 



V2 



1/2 



k 

71 



1 + 



Koje'J 



1/2 



+ 1 



(3.18a) 



Q 



2 

V2 



i\ 2 



_k_ 

71 



(l + tan 2 ^) 1 / 2 -! 
2 \ 1/2 

^ \ V2 



1/2 



l+(^ 



(3.18b) 



The dimensions of a and /3 are m _1 and they are referred to as neper and 
radian, respectively, to indicate their attenuative and phase meanings in wave 
expressions. For lossless media we have a e = 0, a = and the phase constant 
becomes 7 = jj5 = j k. 



2 Recordemos que los valores del factor de atenuacion, tal como se han calculado, vienen 
expresados en nepers/metro y que multiplicados por 8'868 se convierten en dB/m. 
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Equations (3.16) have solutions of the form Ee 1 ^ and He 7 ^ so that the 
instantaneous values for the fields are given by wave equations 

E = ReiEe^-^} = Re{Ee Wu "-^ } = Re{Ee- ai e j{ojt - K) } 

(3.19a) 

H = ReiHe^-^} = RejHe^'^'^} = RejHe^V^*-^} 

(3.19b) 

where the so-called complex propagation vector 7 = 777, (with module 7 and 
direction of the unit vector n normal to the wave-front planes) has been in- 
troduced and r is the position of any point on the wave-front plane so that 
h ■ r = £. 

Equations (3.19) represent waves traveling at a speed given by the phase 
velocity v p 

v v = ^ (3.20) 

which in general, as (3 is given by (3.18a), depends on the frequency (dispersive 
media). 

The penetration factor S is defined as 

5=- (3.21) 

a 

This is the distance at which, due to the attenuation a, the field module de- 
creases from an initial given value to 1/e of this value. 
From (3.7) the following equalities may be deduced 

7-E = (3.22a) 

7 • H = (3.22b) 

7 x E = j> wH (3.22c) 

7 x H = -je ec toE (3.22d) 

From these equations, we see that E, H and n are perpendicular to one 
another and that they form a right-handed system in the order E, H, h. For this 
reason these waves are often referred to as transverse electromagnetic (TEM) 

— * — * 

waves. The magnitudes of E, H are related by 

E 7E „ 
H = — = 4 — 3.23 

Vc J^Vo 

where the quantity i] c , known as the complex characteristic impedance of the 
medium, is given, taking into account (3.13) and (3.17), by 
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Thus, its module and phase is given by 



[l + (^)2]l/4 



(3.25a) 



= tan- 1 | = itan- 1 ^-^ (3.25b) 
Therefore, in general there is a phase shift 9 between E and H. 
3.2.2 Propagation in lossless media 

By particularizing the above expressions for a lossless medium where, e' = s = 
£ r £o, e" = 0, and a = 0, we thus have tanJ^ = 0; 7 = jfc; and 7 = fc = fc n, and 
consequently equations (3.12) and (3.15) simplify to 

V 2 H + fc 2 H = (3.26a) 
V 2 E + k 2 E = (3.26b) 

and the complex characteristic impedance of the medium, (3.25), simplifies to 

= (3.27) 

so that the impedance is real and constant. In particular, when the medium is 
free space, 77 simplifies to the impedance of free space 

\ 1/2 

Mo 



r] = r) = y™j = 120tt (3.28) 

Consequently, in unbounded lossless media, there is no phase shift between E 
and H and the attenuation is null (a = 0). Thus 7 = jk and S = 00 and Eqs 
(3.22) simplify to 

k-E = (3.29a) 

k H = (3.29b) 

k x E = h ujH (3.29c) 

k x H = -jeuoE (3.29d) 



3.2.3 Propagation in good dielectrics or insulators 

In a good dielectric (see Subsection ??) the reactive current predominates on 
the dissipative current and according to (1.91), t&n5d = & e /ue' « 1. In this 
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Figure 3.3i Cuidadojji estan normalizada a 7]q possitive £ traveling fields of a uniform plane 
dissipative medium 




Figure 3.4: Uniform plane wave propagating in the -\~t direction in a lossless medium 
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case, we can develop the complex propagation constant (3.17) to get 



■ ; /, ■, r a/2 • / ,u/2 A jtan<5 d tan 2 5 d 
jk(l - jt<mS d ) 1 = ju(fi s)' [1 + — — 

. , ,,i/2 (-. jta,n5 d \ 

ju (n e ) 1 — (3.30) 



and therefore 



uj{ji Q e') 1/2 tim5 d _ a e /Mo \ 1/2 



2 2 V £ ' 

P ~ fc = w(^ £') 1/2 (3.31b) 



Thus the propagation velocity can be approximated by 

u 1 



" (a^') 1/2 



(3.32) 



From (3.31a) it can be seen that a is small and therefore so is the wave atten- 
uation. Moreover, since a e /ue' << 1, the intrinsic impedance of the medium 
(3.25) is usually simplified to 

V C - V=(^) 1/2 (3.33a) 
9 = (3.33b) 

3.2.4 Propagation in good conductors 

For a good conductor (see Subsection ??) the dissipative current predominates 
on the reactive current and according to (1.94), tanJ^ = a/ws » 1. In this 
case, from (1.93) and (3.13) we have 

1/2 1/2 f ® 

7 = 3 k i l - jtanJd) ~ jfc(-jtan5 d ) = jfc ^— (1 - j) (1 - j)J 
= d + ,)(^) 1/2 " (3.34) 

and consequently from (3.17), 

« = P=(^) 1/2 (3.35) 
Thus the electric field from (3.19a), simplifies to 

E = Re{Ee-^ s e^ t -^ 5 ^} (3.36) 

where 8 

1 / 2 \ 1/2 
5=i=(^f_) (3.37) 
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is the penetration factor (3.21) particularized by a good conductor. Thus, for 
good conductors, the penetration factor <5 has a very low value which decreases 
as the frequency increases. Thus the fields are confined within a very short 
distance from the surface of the conductor. For a perfect conductor, a — > oo 
and <5 = 0. Furthermore, the dielectric constant and the complex impedance 
are reduced to 

Sec = e(l- J —)~ -j- (3.38) 

\ CUE / LO 

and, respectively 




(3.39) 

Thus the phase shift between E and H is 45°. 
3.2.5 Surface resistance 

Let us consider an area element perpendicular to the direction of propagation £. 
Since the wave amplitudes of E and H decrease exponentially according to the 
factor e~ Q £, the complex Poynting vector (1.107), and consequently the mean 
power per unit of area, (1.106), attenuates along the direction of propagation 
by the factor e~ 2a ^. Therefore 

Pav = \ Re{.E x H*} = V av (0)e- 2ai (3.40) 

where V a v{0) is the mean power per unit area at £ = 0. Thus the total power 
per unit area transmitted by the wave to the medium along the distance £ = I 
is given by 

rIP -> 

-- = Pav(0) - V av {l) = P av (0)(l - e- 2al ) (3.41) 
This can be also calculated, according to (1.87), as 

f = ^ = ^d-«- M ) («2> 

This expression for I = oo, or for a distance I such that the magnitude of the 
fields becomes negligible, simplifies to 

^ - ^ \ R ^)El - \ Re{ Vc }H 2 = (3.43) 



since 



(3.44) 



80 



CHAPTER 3. ?? ELECTROMAGNETIC WAVES 



For a good conductor, expression (3.43) simplifies, from (3.39), to 
where R s is the so-called surface resistance 
and 5 is the penetration factor given by (3.37). 

3.3 Group velocity 

So far, we have considered the ideal case of a plane harmonic wave, i.e. one 
in which the wave number and the frequency are fixed. When this type of 
wave propagates through a dispersive medium, the propagation velocity (phase 
velocity) of a harmonic wave depends on its frequency. In practice, the ideal 
situation of a pure harmonic wave which extends to infinity both backward 
and forward in time never arises and, moreover, such a wave could not carry 
information. What in fact happens is that a transmitter emits a given signal 
/(£,£) for a finite period of time that, according to Fourier's theorem, can be 
expanded into a continuous spectrum of amplitudes A^ such that 

/oo 
A^^-^duj (3.47) 
-oo 

When the signal propagates through a dispersive medium, i.e. a medium 
where the phase velocity depends on the frequency, each spectral component 
travels at a different velocity and, as a consequence, the signal will deform as 
it propagates. When, as commonly occurs in practice, the spectrum of the 
signal is narrow 3 and the transmission medium is only slightly dispersive, then 
a single velocity, termed the group velocity, may be assigned to the signal which 
is usually known as a wave group or wave package. The velocity with which 
the envelope or energy of the wave group propagates in the medium is called 
group velocity. To calculate this, let us consider a wave group centered on a 
frequency lu such that A^ ~ except for w = Wo± Aoj/2 (Fig. 7.4). Under 
these conditions, Eq. (3.47) simplifies to 

f(£,t)= I A^-^duj (3.48) 

extended to the values of cj in which A w ^ 0. Given that j3 — f3 (co), it can be 
developed into a Taylor series around the frequency u>q 



/3H = /3M+£ 



{u>-u H) )+— 2 

OJQ 



2o t . ^ 



(3.49) 



WO 



3 Note that a concentration of the field in space docs not imply a concentration in the frequency 
spectrum, hut just the opposite, in accordance with the scale change property of the Fourier transform, 
which indicates that an inverse relation exists between the duration of a signal and its bandwidth. 



3.4. POLARIZATION 



81 



If the dispersive medium is such that the dependence of the phase velocity v p 
on the frequency is so slowly that we can consider (as a good approximation) 
that there exists a linear relation between (3 and w, then (3.49) simplifies to 



P M = 0o + 



dp_ 



(uj - LJ ) 



(3.50) 



where (3 = f3 (w ). 

By substituting (3.50) in (3.48) we get 



H S§L<-oe- 



%i ) [ Aj^-^)^ (3.51) 

J Aw 



which, taking into account (3.48), can be written as a function of f(0,t) in the 
following way 



mt) = f \o,t- 



(3.52) 



This means that, at a point £, the signal has the same amplitude as at the origin 
after a time t = d(3/dcj\ LOg £ and a phase shift given by d(3/du\ LUg Cl> £ — /3 ^. 
Consequently, the velocity at which the signal, and thus its associated energy, 
propagates is 



v„ = 



dt; d(J 

~dt ~ dj3 
d_ 

dp 



U} 



(v p I3) 



dv p 



, a dv P 
1 



d/3/duj 



(3.53) 



If the phase velocity varies slowly with the frequency, then a pulse may 
travel through a dispersive medium a certain distance without a significant 
change. If this condition is not satisfied and the medium is very dispersive 
the shape of signal changes rapidly and the concept of group velocity is not 
longer valid. The sign of dv p /dui determines whether v g is greater or less than 
v p . If the phase velocity v p increases with the frequency, it is termed normal 
dispersion. On the contrary, when v p decreases with the frequency, it is termed 
anomalous dispersion. In an ideal dielectric where v p v p ((3), so that all the 
wavelengths propagate at the same velocity v p = v g , the signal propagates 
without deformation. 



3.4 Polarization 



As the wave equation is a linear differential equation, it fulfils the superposition 
principle and any sum of solutions is also a solution of the differential equation. 
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In particular, let us consider the sum of two plane waves propagating in direction 
z (one with the electric field lying along the x axis and the other along the 
y axis) at identical frequencies but, in general, with different amplitudes (a 
and b) and phases (Si and #2), respectively. Each of these waves, because the 
direction of their electric field does not change with time, is said to be linearly 
polarized, one in the x direction and the other in the y direction. However, 
in an electromagnetic wave the direction of the electric field generally changes 
and traces out an ellipse as the wave propagates 4 . To see this, let us consider 
the total time-varying electric field, which is sum of the two linearly polarized 
waves, given by 

E(z, t) = (ae jSl x + be? s *y)e? { f*- kx '> (3.54) 

Let us determine the time evolution in a plane z = cte of the electric field 
vector resulting from the composition of these two plane waves. We will assume 
a homogeneous, isotropic, lossless medium (although the effects of losses as an 
exponential factor common to all the field components do not influence the 
polarization). 

At the plane z = 0, for example, we have 

E x = acos(ujt + 8i) (3.55a) 
E v = bcos{ujt + 5 2 ) (3.55b) 
E z = (3.55c) 

Using the trigonometric identity for the sum of two angles, solving for cos uit 
and sinoji in terms of a and b, defining 8 = 8\ — 82 as the relative phase 
difference between the two components and after some simplifications based on 
simply trigonometric identities, we find 

E 2 E 2 9F, F, 

— f- + -j7T — -^p 1 cos 5 = sin2 & ( 3 - 56 ) 
a z b z ab 

which is the equation of an ellipse with its major axis tilted depending on the 
value of 8. This means that at a plane z — cte, as the time goes on, the electric 
field delineates an ellipse or, equivalently, that the electric field delineates an 
elliptical helix in the direction of propagation. The resulting polarization is 

referred to as elliptical polarization. The angular velocity of the vector E t = E x 
x + E y y is given by 

ip = "1 = ^an- 1 §0 = E ' Ev -,f vE ' (3.57) 



dt dV E x ' |R 



where ip is 

The sense of rotation together with the direction of propagation define left- 
handed polarized versus right-handed polarized waves, according to the right- 
hand rule: the thumb of the right hand is pointed in the direction of propagation. 



In a unpolarized wave, the vector E is subject to random changes of amplitud and phase 
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Thus, if the fingertips are curling in the direction of the rotation of the electric 
field, the wave is right-handed polarized, and in the contrary case the wave is 
left polarized. 

Particular cases occur depending on the values of a, b, 5, and the polarization 
ellipse may degenerate into a centred ellipse, a circle or a straight line. 

When a ^ b and 5 = mir/2, with m — ±1, ±3, ±5, .. the polarization ellipse 
(3.56) becomes a centred ellipse with the major and minor axis oriented along 
the x,y directions, i.e. 

E 2 El 



— + 



IP 



1 



If a = b, then 



El 



(3.58) 
(3.59) 



which is the equation of a circumference. 

When 5 = ±mir, with m being an integer, the equation (3.56) becomes 



a b 







which represents the equation of a straight line 



Ey — ^f — E x 

a 



(3.60) 



(3.61) 



intersecting the origin. The wave is then linearly polarized and the components 
of E are 



E x = acos(uit~kz) 

E y = bcos(u>t — kz ± m7r) 

The angle of the slope with the x axis is 



tan ip = tan 



E v 
E x 



(-1)' 



(3.62a) 
(3.62b) 



(3.63) 
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Chapter 4 



Reflection and refraction of 
plane waves 



In the previous chapter, we studied the characteristics of harmonic plane waves, 
and now consider what happens when such waves reach the interface (assumed 
to be plane and indefinite) separating two linear, nonmagnetic, homogeneous 
and isotropic dielectrics having different electromagnetic characteristics. The 
change in the constitutive parameters, as the wave passes from one medium to 
the other, is assumed to take place in an electrically very narrow region with 
a thickness much less than A. In general, when a wave propagating through 
a medium strikes the interface (incident wave), part of its energy is reflected 
and propagates through the same medium (reflected wave), while another part 
is transmitted to the second medium (transmitted, or refracted wave). The 
characteristics of reflected and transmitted waves can be calculated from those of 
the incident wave by forcing the total field on the interface to fulfil the boundary 
conditions. We will consider first the simplest case of normal incidence, i.e. 
when the interface is perpendicular to the propagation direction of the wave, 1 . 
and then the more general case of oblique incidence. This study has extensive 
applications in optics where the interface of many optical devices, such as lenses 
and fiber-optic transmission lines, has a radius of curvature much larger than 
the wavelength of the incident wave. Thus the interface can be considered quite 
accurately as a plane interface. In the following, with no loss of generality, we 
will assume the interface to be parallel to the xy plane. 



La incidcncia normal tiene muchas analogfas con la h'neas de transmision que se estudiaran en 
el capftulo Tal 
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Figure 4.1: Poner los vectores de pynting V El subindice de campo electrico 
incidente ponerlo mejor 

4.1 Normal incidence. 

4.1.1 General case: interface between two lossy media 

normally incident from a lossy media, characterized by the parameters fi ,£i = 
e\ — je'l, u\ to the surface of another one with different constitutive parameters 
Mo> £ 2 = 4 -je'2,a 2 

Considering two semi-indefinite lossy media that are separated by the plane 
z = 0, see figure 4-1, let us assume that a harmonic plane wave propagates 
through the first medium in the positive sense of the z axis with the electric field 
parallel to the x axis. The wave impinges with normal incidence on this plane. 
Due to the discontinuity of the constitutive parameters, u, = u, , s C i = e\ —je", 
and c>i where subindex i (i — 1,2) refers to medium 1 or 2, part of the wave 
is propagated through medium 2 and part is reflected back through medium 
1. Therefore, the total field in medium 1 (where z < 0) and medium 2 (where 
z > 0) is given by 

Medium 1 

E xl = El ie - aiZ e- J ^ z +E r xl e aiZ e 3 ^ z =E l xl e-~^ z +E T xl e~<i z 

(4.1a) 

pii jpr pii pir 

Vcl Vcl Vcl Vol 

(4.1b) 

Medium 2 

E x2 = E x2 e- a - z e-i^ z =El 2 e-^ z (4.1c) 

Vc2 Vc2 
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The ^^superindices?? %. 7", and t indicate the incident wave {medium 1), the reflected wave 
(medium 1) and the transmitted wave (medium 2), respectively. TllG millUS Sign 101 the 

reflected wave of the magnetic field is associated with the fact that the Poynting 
vector of the reflected wave propagates in the — z direction. In these expressions, 
Vd = V Mo/ £ eci represents the impedance (3.24) of medium i while j i is the 
complex propagation factor (3.17), 



7; 



■3 Pi 



(4.2) 



where on and j5 i are the attenuation and propagation constants (3.18a)and 
(3.18b), respectively 



Pi 



a; 



V2 ' 



+ i 



1 + (<7 e /<) - 1 



1/2 



1/2 



(4.3) 



(4.4) 



The time dependence of the fields is achieved by adding the factor e JW * to (4.1). 

For each instant of time, by imposing the boundary conditions in the plane Z — (2.182b) 
onto the tangential components of E and H , 



we obtain 



to 

tj x2 



H 



En — E<n 
Hu — H 2 t 



uj = uj = UJ 

t lE 1 x1 = (1 + T L )E l xl 

Vcl, 



Vc2 



where is the reflection coefficient in the plane z = defined by 

r L = Vc \ Vcl =\r L \ e ^ 

Vc2 + Vcl 

and Tl is the transmission coefficient in the same plane, defined by 

2^c2 



Vc2 + Vcl 

1 + T L = \T L \eP Vl 



(4.5a) 
(4.5b) 



(4.6a) 
(4.6b) 
(4.6c) 
(4.6d) 

(4.6e) 



(4.7) 

(4.8a) 
(4.8b) 



If there is an impedance adaptation (f] c2 = f] cl ) then there is no reflected 
wave, and so all the incident energy is absorbed by the second medium. 



88 CHAPTER 4. REFLECTION AND REFRACTION OF PLANE WAVES 



From (4.1a) and (4.6b), the total electric field in the first medium can be 
expressed as 

E xl = E l xl e- aiZ e-^i z (l + r L e 2aiZ e 2 i^ z ) (4.9) 
= E xl e- aiZ e-i^ z (l + T{z))=E xl e-~<i z {l + T{z)) (4.10) 

where T{z), defined as 

T(z) = r L e 2aiZ e 2 ^ z (4.11) 

is the reflection coefficient in the plane z = z. Similarly, for the magnetic field, 
we have 

H yl = ^e-^ z (l-T(z)) (4.12) 

The impedance associated with the total field at a coordinate point z in the 
first medium is defined as 



{ \ — Exl 

Hyl 



1 + V c2 -Vci tanh( 7l z) 

Vcl l^fjz) ^ cl -, c2 tanh( 7l z) (413) 



The impedance T] inp (z) is continuous through the interface, because the tan- 
gential components E x \ and H yl are similarly continuous, while the reflection 
coefficient T is discontinuous. 

4.1.2 Perfect/Lossy dielectric interface 

In the particular case in which the first dielectric is perfect, i.e. lossless (ui = 
and e" = 0, z\ = s' 1 , / y 1 = jk\), the characteristic impedances reduce to 

m = M (4-i4) 
V £ i 

and the coefficient of reflection (4.7) at the interface (z = 0) becomes 

\ — j £ ec2 

r L = — (4.15) 

V 61 

Since o\ = and e" = 0, it follows that a.\ — and therefore 

E xl = Ei ie - jklZ (l + T L e 2 ^ z ) = Ei ie - jklZ (l + T(z)) (4.16) 

and 

H yl = ^e-^ z (l-T(z)) (4.17) 
m 

where 

Y{z) = T L e 2 ^ z (4.18) 
The input impedance (4.13) simplifies to 

lmpK 11 m -?7 c2 tan(fc 1 z) 
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4.1.3 Perfect dielectric/Perfect conductor interface 

Another particular case arises when the second medium is a perfect conductor 
(t] 2 = 0) and therefore Tl = and Tl = —1. Then the helds in the first medium 
are 



E x i — 

Hyl = 



K 



(1 - e 2jfclZ ) = E xl (e~ jklZ - e 3klZ ) = -2jE xl sin((42B)i) 



+ e> 



2^icos(fciz) 



(4.20b) 



4.1.4 Standing waves 

It is well known that two waves with the same frequency that are propagating in 
opposite directions interfere and form what are termed standing (or stationary) 
waves. To examine this concept, let us first consider the case in which the 
first medium is lossless, and then analyse the case in which the first medium is 
dissipative. 



a) Lossless case 

For the first medium, the expression of the total electric field is, 

E x i = E xl e~ jklZ + E r xl e jklZ = (1 + T L )E xl e- jklZ + T L E l xl {e^ z - e~ jklZ ) 
= T L E xl e-^ z + \T L \ Bj. 1 (e , ' ( * 1+fciz) - e j ^ L - klz) ) 
= \T L \ Ei ie j ^ L - kl ^ + 2 \T L \E xl sin(fci«) e ^ L+ * /2) (4.21) 

By including the time dependence, and assuming an initial phase tp = 0, we 
obtain the following expression for the total field 

E xl (z,t) = \T L \E l 0xl cos(wt - ktz + * L ) - 2 \T L \E l 0xl sm(k lZ ) sin(wi + § L ) 

(4.22) 

where the first summand of the second member corresponds to a wave that is 
propagating, while the second summand represents a standing wave, i.e., one in 
which the mean energy transported by the wave is null. The amplitude of the 
propagating wave is determined by the coefficient of transmission, while that of 
the standing wave depends on the coefficient of reflection. The envelope of the 
equation (4.22) is termed the diagram of the standing wave. If the coefficient 
of transmission T is null (which occurs when the second medium is a perfect 
conductor) the wave of the first medium becomes a pure standing wave. 
From (4.16) and (4.17) the magnitudes of the fields are 



E 0xl = E* 0xl \l + T(z)\ 

H yl = ~ 



E 0xl I 1 



T(z) 



(4.23a) 
(4.23b) 
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The maximum values of Eq x \ (the minima of Hq v1 ) are given by 

E 0xl (z) inax = E l 0xl +E r 0xl (4.24) 



at the coordinate points 

<f>L + 2mr 

^max ^ 



n = Q,l,... (4.25) 



and the minimum values of Eq x (the maxima of Ho y ), assuming Eq x1 > Eq x1 , 
are given by 

E 0x {z) min = Eq x1 - Eq x1 (4.26) 



at the points 

$ L + (2n + l)7T 



2ki 



; n = 0,l,... (4.27) 



Ratio of the standing wave 



The relation between the maximum and minimum values of the diagram of the 
standing wave is called the ratio of the standing wave, and is described by 

SWR = ^Oxl(z)max = goal + E Qxl = 1 + l r ( Z )l = 1 + \ T ' l\ ^ ^ 

E 0x i(zUn Ei xl -E r 0xl l-\T(z)\ 1-\T L \ 
Its value ranges from 1 (no reflected wave) to infinity (pure standing wave), i.e. 

1 < SWR < oo (4.29) 

b) Lossy case 

In this case, the expression of the total electric field in the first medium is 

E xl (z) = E l xl e-^ z e-3^ z + E r xl e aiZ e j ^ z = E xl (e~^ z + IW 1 *) 

= T L E xl e-^ z +T L E xl {e'^ z - e^ 2 ) (4.30) 

and, by including the time dependence, we get 

E xl (z,t) = \T L \E z 0xl e- a ^cos(ujt- p lZ + * L ) + 

2 \T L \ E l 0xl sinh( 7l z) cos(wt + $ L ) (4.31) 

In these media, it makes no sense to define the SWR parameter because the 
maxima and minima are not constant. 
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4.1.5 Measures of impedances 

Assuming that the first medium is lossless, from (4.27) the first field minimum 
occurs at </> L +2fc 1 z min = ir. Consequently, we can determine the phase angle <f> L , 
assuming that k\ is known and that z m i n is determined experimentally (by using 
an appropriate device to detect the first field minimum). If fciis not known, it 
can be calculated from the distance between two consecutive minima. 

The value of \Tl\ can be found from the ratio between the maximum and 
minimum field values E 0xlmax /E 0x i min = SWR = (1 + |r L |)/(l - |r L |). 

Note that if the incident wave has an amplitude of one, then \T L \ is identical 
to Eq and thus we need to measure only this amplitude. Thus r\ 2 is determined 
from this information and from expression (4.7). 



4.2 Multilayer structures 

Let us now consider the normal incidence of an electromagnetic wave on a 
structure in which there are more than two media separated by parallel planes. 
To simplify the analysis, we consider the case of three lossless dielectrics, as 
shown in Fig. tal. The generalization to more media, including the possibility 
of losses, is straightforward. Clearly, for a wave that is propagating to the right 
in medium 2, the problem is analogous to the two-layer cases discussed above. 
Therefore, the coefficient of reflection in the z = plane is 

T 23 = = T(z = 0) (4.32) 

V3 + V2 

where subindex 23 refers to the surface that separates medium 2 from medium 
3. Particularizing (4.13) for z = — I we have the load impedance rj L 

VL = Vinpiz = -0 = % 1 _ T23e -2 jk l ( 4 ' 33 ) 

Taking into account that i] inp (z) is continuous at an interface, the coefficient 
of reflection (4.7) at z = —I, becomes 

T L =T(z = -l) = ^-^ 
Vl+Vi 

by introducing (4.33) into this equation and then operating, we get 

r 12 + r 23 e- 2 ^ 



1 + r 12 r 23 e- 2 ^ 

where 



(4.34) 



r« = (4-35) 
Vj + Vi 

Thus, for an electromagnetic wave with an amplitude of one, impinging 
normally from the first medium onto the structure, the amplitude of the reflected 
wave is given by (4.34) 
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Quarter-wave layer 

For a quarter-wave layer, I = A/4 [e~ 2: > kl = —1), equation (4.34) becomes 



i-r 12 r 



12 1 23 



(4.36) 



Thus, to transmit the incident energy completely (adaptation of impedance), 
the coefficient of reflection must be null, and so 

Ti 2 - T 23 = (4.37) 

Taking into account equation (4.35) we have 

f)i = Vmm ( 4 - 38 ) 

as a condition for impedance adaptation to exist. 
Half- wave layer 

For a half-wave layer, i.e. I = A/2 (er 2 ^ kl = 1), expression (4.34) is reduced to 

For impedance adaptation to exist, the following must be fulfilled 

Ti2 + r 23 = (4.40) 
By replacing the coefficients by the values given in (4.35), we have 

(% = m) (4.4i) 

and thus = irrespectively of rj 2 . Thus any material with a thickness of A/2 
is adapted so long as the impedances of media 1 and 3 are the same. 

4.2.1 Stationary and transitory regimes 

The above analyses are valid for monochromatic waves in a stationary regime. It 
should be noted that such a regime is the limit of a transitory process involving 
multiple reflected and transmitted waves within media 1 and 2. To illustrate 
this limit process, let us consider the normal incidence of a wave that impinges 
upon a structure formed of three perfect dielectrics, as shown in Fig. tal. From 
the process of multiple reflections and transmissions, we find that in medium 1 
a reflected field is given by 

Ki = 3§xi(ri2 + T 12 T 23 T 21 e- 2 ^ d + T 12 Y 2 2 . i T 2l T 21 e- A ^ d 

+T 12 Tl 3 T 2 21 T 21 e- G ^ d + ...) (4.42) 
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Figure 4.2: Poner sistemas de ejes 



Observing the second member, we can see that the summands following 
the first one constitute a geometric progression of common ratio r 2:L r 2 3e~ 2: ' fe2 ' i . 
Thus the coefficient of reflection can be written as 



r L = r 12 + 



Ti 2 r 2 3T 2 ie 



-2jk 2 



i - r 23 r 21 e- 2 ^ 

which, taking into account the equalities 

r 2 i = — r i2 
r 12 = i + r 12 
T21 = 1 — T12 

is reduced to equation (4.34). 



(4.43) 

(4.44) 

(4.45) 
(4.46) 



4.3 Oblique incidence 

As a more general case than the normal incidence, let us now consider the 
oblique incidence of a plane wave on a plane interface separating two media. In 
general, in medium 1 there exists an incident and a reflected wave, while the 
transmitted (also called refracted) wave is in medium 2. To study the oblique 
incidence we will use the geometry shown in Fig.??, where the waves have been 
represented, as usual, by arrows (called rays) in the direction of propagation. 
These rays are perpendicular to the equiphase planes (wavefronts). The oblique 
incident has extensive applications in optics where the interface of many optical 
devices, such as lenses and fiber optic waveguides, has a radius of curvature 
much larger than the wavelength of the incident wave. Thus the interface can 
be considered very approximately as a plane interface. 

In principle, we make no assumption that the three rays are coplanar, al- 
though they are shown as such in Fig Tal. The plane of incidence is defined 
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by vector 7 1 and by the z axis. Let us assume that 7* is in the plane y = 
and forms an angle $i with the z axis. In the general case of two lossy media, 
the electric fields of the incident, reflected, and refracted waves can be written, 
respectively, as 

E l = Re{^e (j ' wi '-^} (4.47a) 
E r = Re{i%e^ r *-^} (4.47b) 
& = Re^e**"'* - ^} (4.47c) 

In z = 0, the tangential component of the electric field must be continuous, 
and thus we have 

K+E:=El (4.48) 

so that 

Re^e^*"^} + Re^e^*-^} = Re^e^'*"^} (4.49) 

A similar relation must be fulfilled between the components of the fields with 
respect to the y axis. These conditions can be satisfied only if 

J = cj* = u r = to (4.50) 

and 

7* • f '= 7 r • f '— 7 • r (4-51) 

Since 7* lies on the y — plane, from (4.51) it follows that 7^ = 7^ = 0, 
signifying that the reflected and refracted waves are coplanar with the incident 
wave. Thus we have 

[xsin^i + zcos#;] (4.52a) 
[a; sin 6 r + z cos 8 r ] (4.52b) 
[x sin 8 t + z cos OA (4.52c) 

c 

where JV C is the complex index of refraction of the medium, such that 

7 = ^ (4.53) 

c 

By substituting (4.52) in (4.49), and by making the coefficients of x equal, 
we get 

6i = 8 r = 6 (4.54a) 
N cl smd = N c2 sin8 t (4.54b) 

These equations, together with the coplanarity of the rays, constitute Snell's 
laws. 



f 


• r = 


N cl LU 

c 


-*r 
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■ r = 


N c1 uj 
c 


f 


■ r = 


N c2 LU 



4.4. INCIDENT WAVE WITH THE ELECTRIC FIELD CONTAINED IN THE PLANE OF INCIDENCE95 



For lossless media Eq. (4.54b) simplifies to 

N lS in9 = N 2 sin0 t (4.55) 

where 

Ni = — = (At ri e ri )' (4.56) 

or, for nonmagnetic media, 

Ni = (e ri )^ (4.57) 

Next we study the relations between the amplitudes of the incident, trans- 
mitted and reflected waves by making use of the boundary conditions at the 
interface between the two media. For this we will assume lossless media al- 
though the generalization to lossy media is straightforward 2 . Let us analyze 

the problem in two stages, firstly where the electric field E oscillates in the 
incidence plane, and then where it oscillates perpendicularly to the same plane. 
Any other case can be considered a superposition of these two situations. 



4.4 Incident wave with the electric field con- 
tained in the plane of incidence 

From the continuity of the tangential components of E and H (Eqs. (4.5a) and 
(4.5b), we obtain (Fig. 9.2) 

Ef|Cos6» + E||cos6»-E||Cos6» t = (4.58a) 

— {E\-E\)-—E\ = (4.58b) 

Vi " 11 V2 

where the subindex || indicates that the physical magnitude in question lies in 
the incidence plane and 

E\ = & m e-&-f (4.59a) 
E\ = ELe-^- (4.59b) 



From (4.58) we find that 



" El V2 cos "t + Vi cos " 



-Ey 2ri 2 cos9 



El r\2 cos ®t + ?7i cos 9 



(4.60b) 



Where Ty and ry are the coefficients of reflection and transmisison, respectively. 
If medium 2 is a perfect conductor (rj 2 = 0) then T\\ = —1. 
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Figure 4.3: cuidado con superindices y subindices... 

Li>lF or lossless non-magnetic materials ?? (}J,^ — — I^q) such that 

r] 2 v 2 Ni sin6> t 
JV12 = — = = -rr- = . - 
Vi v i N 2 sm9 



cuidado notacion dc vol 
fasc. 



(4.61) 



with N\2 being the ratio of the indices of refraction of medium 1 and medium 2. expressions 

(4.60a) and (4.60b) are reduced to 



EJ 

Em 



tan (9 t - 9) 
tan (9 t + 9) 

2 sin 9 t cos 9 



E\ sin {fit + 9) cos (9 t - 9) 



The total electric field Em in medium 1 is given by 



where 



^1 -i -r 

Eu = E i + Ei 



Ey = E\cos9x - E\sin9z 
En = Eli cos 9x + Eii sin 9z 



(4.62a) 
(4.62b) 

(4.63) 

(4.64a) 
(4.64b) 



Thus we have 



E J = cos (Ej,| e-^'-*' + e^'^) 4 + sin (e^, e"^^ - E*,, e^' 7 



z 

(4.65) 



If the medium is lossy, we must replace {r},jk} by (?7 C ,7) 
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that is, 



E\ x = E^cosO (e-^- p + T^~ kT -r) (4.66a) 
E\ z = E^sme^e-^-r-e-^ (4.66b) 

Taking into account that 

k r ■ f = — k r z cos 8 + k r xs'm6 (4.67a) 
k l -r = tfzcosd + k l xsm6 (4.67b) 

and by substituting these equations in (4.66), we find that 

E\ x = El n er^ ixs ' ln6+zcose) {l + T\\e 2 ^ zcosd ^cos6 (4.68a) 
E\ z = -El n e-' jk ' {xs{n6+zcos ^ (l _r ||e 2jfe^co S e^ gin6l ( 468b ) 

When the time factor e Ja; * is introduced, the term e^"* - -'* xslnd ) represents 
a wave that is propagating in the direction of the x axis, while the term 

^-jk'z COS0 +T ^ e jk i z cosO^ e jut ^gg) 

llll llll llF rom (4.68a) or the corresponding one from (4.68b) gives us the 
superposition of two waves that are propagating with respect to the z axis, but 
in opposite directions. In other words, a stationary wave overlies a traveling 
one such that the energy that is transported in direction z, from medium 1 to 
medium 2, is transported by the traveling wave.????. For the case of a perfect 
conductor, = — 1 and there exits only a standing wave along the z axis. 

The total magnetic field H ± in medium 1 is 



H\ = Hi ± e~ jki - f + H5 ± e- iUr - f = ^ e'^ - E^ e - J '* v ) y = 

Jji ^ e ~jk\x sin 6+ z cos 6) ^ _ e 2 3 k ' z cos (4 70) 

where the symbol _L indicates that the magnitude in question is perpendicular 
to the incidence plane 

In the case of a perfect conductor, it is straightforward to show that there is 
no energy flow towards z, but there there is towards x, as the mean time value 
of Poynting's vector towards z is zero. 
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Figure 4.4: cuidado el reflejado tiene mal el sentido del campo electrico 

4.5 Wave incident with the electric field perpen- 
dicular to the plane of incidence 

As above, the continuity equations are used for the tangential components of E 
and H and thus (Fig.9.3): 

E ± +E r ± = E ± (4.71) 

-*i -*r -*t 

E± ~ E± cos9 = ^cos0 t (4.72) 
Vi f)i 

where the subindex _L indicates that the physical magnitude in question corre- 
sponds to the case in which the electric field of the incident wave is perpendicular 
to the plane of incidence. 

By resolving the two Eqs. (4.70) and (4.71), we get 

= E 1= y 2 co S 9-r ll co S t 
E\ t) 2 cost) + r\ x cos 6*4 

E z j_ T)2 cost) + r\ x cos 6*4 

where the parameters T± = El/Eland t± = E ± /E l ± are the coefficients of 
reflection and transmission, respectively. 

For lossless non-magnetic materials, Eqs. (4.73) are transformed into 

El Binfr-fl) 

El sm(e t + 9) (4 ' /4aj 

El 2sing t cosfl 

El sin(0 t +0) 1 J 
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By operating in a similar way to that described for the case of El , we arrive 
at the following for the total electric and magnetic fields in a lossy medium 1 

E\ = E l 0± e- jkixsine (e~^ izcose +r ± e jkZzcos0 ^)y (4.75a) 
Hj = cos<9 (eT&- f - T±e- iUr - F ^ x 

sin0 (e- j *- ? + T ± e- iir - f } z (4.75b) 

As in the case of E\\, the field behaves as a travelling wave towards x and 
as a travelling wave overlying a standing one towards z. The formulas (4.62) 
and (4.74) are known as Fresnel's formulas, which give the relations between 
the amplitudes and phase of the incident, reflected, and tranmitted waves. 
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Chapter 5 



Electromagnetic 
wave-guiding structures: 
Waveguides and 
transmission lines 

5.1 Introduction 

There are many engineering applications in which it is necessary to use devices 
to confine the propagation of the electromagnetic waves in order to transmit 
electromagnetic energy from one point to another with a minimum of interfer- 
ence, radiation, and heat losses. Although such transmission systems can take 
many different forms, a common characteristic is that they are uniform. That 
is, their cross-sectional geometry and constitutive parameters do not change in 
the direction of the wave propagation z for wavelengths numerous enough to 
make border effects negligible. In general, any device used to transmit con- 
fined electromagnetic waves can be considered a waveguide; however, when the 
transmission device contains two or more separate conductors the term "trans- 
mission line" is generally used instead of "waveguide". Figure (5.1) shows the 
cross-sectional shape of some guiding transmission systems: two-wire trans- 
mission lines; coaxial transmission lines formed by two concentric conductors 
separated by a dielectric; two hollow (or dielectric-filled) metal tubes of rectan- 
gular and circular cross section (i.e. a rectangular and a circular waveguide); 
two planar transmission lines (the stripline and microstrip); and two dielectric 
(without conducting parts) waveguides: the circular dielectric waveguide (or 
homogeneous dielectric rod) and the optical fiber. 

In hollow conducting pipes waves propagate within the tube, whereas in 
transmission lines formed by two or more conductors, the waves propagate 
in the dielectric medium between the conductors. In homogeneous dielectric 
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waveguides the Geld decays exponentially away from the dielectric in the trans- 
verse plane, and consequently the electromagnetic waves are conhned mainly 
within the dielectric medium. Optical fibers, used mostly at optical wavelengths, 
consist of a cylindrical core surrounded by a cladding and are usually circular 
in cross section. The light is essentially confined to the core (which has a larger 
refractive index than the cladding) by total internal reflection as it propagates 
along the fiber and the wave is confined without need of any conducting walls. 

The choice of a specific transmission system depends on the application 
and should take into account aspects such as frequency range, losses, power- 
transmission capacity, and production costs. For example, the two-wire trans- 
mission lines, which are usually covered by polyethylene, are relatively inexpen- 
sive to manufacture, but radiation losses (mainly at discontinuities and bends) 
make them inefficient for transferring electromagnetic energy farther than the 
lower range of microwaves. Coaxial lines and hollow metal pipe waveguides are 
more efficient than two-wire lines for transferring electromagnetic energy be- 
cause the fields are completely confined by the conductors. For the transmission 
of large amounts of power at high frequencies, waveguides are the most appro- 
priate means. In a coaxial cable, significant wave attenuation occurs at high 
frequencies because of the large current densities carried by the central conduc- 
tor, which has a relatively small surface area. On the other hand, waveguides are 
intrinsically dispersive and consequently incapable of transmitting large band- 
width signals without distortion. However, coaxial lines can guide signals of 
much higher bandwidths than waveguides can. 

As shown in the next chapter, the dimension of the cross section of a 
waveguide is related to the wavelength of the guided wave. Thus, for very low- 
frequency waveguides the cross section would be too large and thus impractical 
for frequencies lower than 1 GHz. On the other hand, at optical frequencies 
the size of a metal waveguide must be too small (in the range of the urn) and, 
moreover, at these frequencies the study of the interaction of the electromagnetic 
field with the metal walls requires of quantum mechanical theory. 

As a result of the development in solid-state microwave and millimeter tech- 
nology, planar transmission lines are used instead of waveguides in many ap- 
plications because these lines are inexpensive, compact, and simple to match 
solid-state devices using printed-circuit technology. Planar lines allow different 
configurations, usually including a dielectric substrate material with a ground 
plane and one or more conducting strips on the upper surface. The most com- 
monly used of these are striplines and microstrips, which are briefly described 
in Section ??. 

The field configurations that can be supported for any guiding structure 
must satisfy Maxwell's equations and the corresponding boundary conditions. 
The different field distributions that satisfy this requirement are termed modes. 
Although the electromagnetic field distribution in ideal guiding transmission 
systems (composed of perfect conductors separated by a lossless dielectric). can 
be expressed as a superposition of plane waves, the study of the propagation is 
greatly simplified when we seek other kinds of solutions called transverse mag- 
netic (TM) modes, transverse electric (TE) modes, or transverse electromagnetic 



5.2. GENERAL RELATIONS BETWEEN FIELD COMPONENTS 



103 




(5) (6) (7) (8) 



Figure 5.1: Examples of waveguides and transmission systems: (1) Two wire 
transmission line (2) Coaxial transmission line (3) Rectangular waveguide (4) 
Circular waveguide (5) Stripline (6) Microstrip (7) Circular dielectric waveguide 
(8) Optical fiver cable. 

(TEM) modes. These terms indicate that, in the direction of propagation, the 
TM modes have no magnetic field component, the TE modes have no electric 
field component, and the TEM modes have neither electric nor magnetic field 
components. In practice, these modes form a complete set of orthogonal func- 
tions and, hence, any propagating electromagnetic field in the guiding structure 
can be expressed as a linear combination of these modes. As discussed below, 
there are two important properties that distinguish TEM from TE and TM 
modes: 

1) TM and TE modes have a cutoff frequency below which they cannot 
propagate, which depends on the cross-sectional dimension of the guiding struc- 
ture. 

2) TEM modes cannot exist within a waveguide formed a single perfect 
conducting pipe while transmission lines can in general support TE, TM and 
TEM modes. 

In this chapter, we present some general aspects of the propagation of time- 
harmonic electromagnetic waves in guiding systems formed by perfect con- 
ductors and only one homogeneous lossless dielectric in which the guided field 
propagates. Nevertheless, the results can serve as a basis for structures in which 
the cross-section contains more than one dielectric medium. The effect of lossy 
media is analyzed in the final section. The study of some specific geometries is 
left for the next chapter. 

5.2 General relations between field components 

Let us assume that a time-harmonic wave propagates along the z-axis, in the 
+z direction, in a lossless guiding transmission system. Thus the dependence 
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on z and time t is given by the factor e^^^and the fields are of the general 
form 

E oe j{u>t-f3 s z) ] ^ (- g e j ut 



where E = E e _: " 3 9 z and H = Hoe~^s z and /3 is the wavenumber of the 
guided wave. Because the_geometry and constitutive parameters do not change 
along the z-axis, E and H are functions only of the transverse coordinates. 

To determine E and H, we will first show that it is possible to express 
their transverse components, E t and H tl in terms of their z-components, E z 
and H z . For this, we divide the three dimensional Laplacian operator V in the 
homogeneous Helmholtz wave equations (3.26) into two parts. One part, d 2 /dz 2 , 
acts only on the axial coordinate, z, and the other, V 2 , on the transverse ones 
only 1 , i.e. 

V2 = £ + V * (5 ' 77) 
Since d/dz = —j(3 g , the wave equations can be written as 

(V 2 + ^ 2 ){|} = (V? + / 1 2 ){|}=0 (5.78) 



where 



h 2 = k 2 -/3 2 g (5.79) 



and k — ui(/ie)2 is the wavenumber for the wave propagating in an unbounded 
medium of the matter which fills the transmission system. By particularizing 
(5.78) for the z field component, we have 

(v 2 + /> 2 ){h}=° ( 5 - 8 °) 

This equation, when solved together with the boundary conditions of a given 
structure, has solutions for an infinite but discrete number, to, of characteristic 
values (eigenvalues) h m , i.e. 

(V 2 + /4 l ){^ m )=0 (5.81) 

where 

h 2 m = k 2 - f3 2 gm (5.82) 

with E zm , or H zm being the corresponding functions characteristic (eigen- 
functions) which satisfy the equations (5.81) and the corresponding boundary 
conditions, which are determined by the geometry of the system. 



^For example in Cartesian eoordinatcs we have V = Vt + zJ^ where Vt = "t~ ~§y^ so 

tl-t V? = + #j. 



5.2. GENERAL RELATIONS BETWEEN FIELD COMPONENTS 



105 



Now we are going to demonstrate that, once equation (5.80) has been solved, 
we can obtain E t or H t from E z and H z . From Maxwell's equations (1.67c) 
and (1.67d), in a sourceless region, we have 

VxE = -jujuH (5.83a) 
V x H = jueE (5.83b) 

The transverse components of these equations can be written as 

(V x E) t = V t xE z + V z xE t = -ju^Ht (5.84a) 
(V x H) t = V(XH z + V z xH t = jujeE t (5.84b) 

Thus, as E z and H z are assumed to be known, we have a system of two equa- 
tions and two unknowns, E t and H t , the solutions to which are 



H 



= fa (cjs\7 t xE z - /3 g V t H z ) 



E 



(5.85a) 
(5.85b) 



According to (5.85), once the z components of the fields are known, the trans- 
verse components can also be calculated. Moreover, in ideal guiding ^struc- 
tures?? , we can express any field propagating in the homogeneous guiding 
transmission structure as a linear superposition of TE, TM and TEM waves or 
modes. Clearly, it is not possible to find specific expressions for the field dis- 
tribution of any of these modes without previously knowing the geometry and 
characteristics of the transmission system. However, as shown below, we can 
study some of their general characteristics. 

5.2.1 Transverse magnetic (TM) modes 

Let us first consider TM modes so that H z — in (5.85). Thus we have 



E t = - J i£V t E z = V t i 



_ V7 18E,_ 



dz 



t^TM 



(5.86a) 



Ht = xE z = -iffz x V t E z = fixE t = x E t 



(5.86b) 



where we have defined the scalar potential for the TM modes, $tm, as 



TM = pr- 



1 dE z 



(5.87) 



To obtain (5.86b), we have used the equality V f xE z = -zx V t E z and defined 
the frequency-dependent quantity Ztm, as 



Z TM = dr = T 1 



(5.88) 
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where r/ = (/i/e) 5 is the intrinsic impedance of the dielectric that fills the trans- 
mission system. The quantity Ztm, which has the dimensions of impedance, is 
called the wave impedance for the TM modes. From (5.86b), we can see that 
E t , H t , and z form a right-handed system when the wave propagates in the 
z-positive direction. 

Thus, from (5.86a), in TM modes, E t can be written as the gradient of a 
scalar function $™ ■ This result could have been obtained by simple reasoning 
from Faraday's law (5.83a), taking into account that, since H has only trans- 
verse components the same is true for V x E. Therefore, from Stokes' theorem, 
we have 



(V x E) • zds = d> E-dl = 
s Jr 



(5.89) 



where S is a transverse surface normal to the z axis. But E z cannot contribute 
to the line integral because the integration path T lies on the transverse plane. 
Therefore 



E-dl 



E t -dl = 



(5.90) 



which implies that E t is conservative and, therefore, can be written as the 
gradient of a scalar function 3>tm- 



5.2.2 Transverse electric (TE) modes 

For TE modes, from equations (5.85), with E z = 0, we have 



y{2 V t-Tl 2 



Y7 1 dH z 

v */i 2 dz 



t^TE 



(5.91a) 



E t = -W^t xH z = ^zx V t H z = -fz xH t = -Z TE z x H 



where 



1 dH z 



TE - TPST 



is the scalar potential for TE waves and 



ry UJtl Tlk 

ATE - 7T - 7T 



(5.91b) 
(5.92) 

(5.93) 



is the wave impedance for the TE mode. From (5.91b) we can see that E t , H t , 
and z form a right-handed system when the wave propagates in the z-positive 
direction. 

The fact that, according to (5.91a), H t can be expressed as the gradient of 
the scalar function $te can be explained by Ampere's law, (5.83b), following a 
reasoning similar to that used in the case of TE modes. 
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5.2.3 Transverse electromagnetic (TEM) modes 

For TEM modes, since E z — and H z = 0, substituting these values in (5.85), 
we can get no null or trivial solutions only if h = 0. Consequently, from (5.79), 
for TEM modes, we have 

2 g = k 2 = uj 2 fj,e (5.94) 

This means that a TEM mode in a transmission system has the same propa- 
gation constant as a uniform plane wave traveling in the unbounded dielectric 
between the conductors. Since h — and E — E t and H — H t , (5.78) reduces 
to 

V 2 E = V 2 E t = (5.95a) 
V 2 H = W 2 H t = (5.95b) 

Thus the distribution of the electric and magnetic fields on a transverse 
plane satisfies the same bidimensional Laplace's equation as for the static fields. 
This means that, for TEM modes, on a transverse plane, E is conservative, and 
derivable from a scalar function # by means of the gradient function, i.e. 

E = V* (5.96) 

Hence, the electric field distribution in the cross-sectional plane has the same 
spatial dependence as the electrostatic field created by static charges located on 
the conductors of the transmission system. Consequently, a TEM mode cannot 
exist within a waveguide formed by a single perfect conducting tube of any 
cross section since no electrostatic field can exist within a sourcesless region 
completely enclosed by a conductor. When two or more separated conductors 
exist, as for example in coaxial, two-wire or stripline transmission lines, TEM 
waves can be propagated along the dielectric separating the conductors. 
It is straightforward from (5.84b) that 

Et = -Ztem zxH t = -nzx H t (5.97) 

where Ztem 

i 

Ztem = ti={^Y (5.98) 

is the wave impedance for the TEM mode, which coincides with the character- 
istic impedance r\ of the dielectric that fills the transmission system. cuidado en lo de lineas: usar 

Now we will demonstrate that, for TEM modes, Maxwell's equations can be o no negritas..?- 
used to derive a pair of coupled differential equations which enable us to study 
the propagation of these modes in transmission lines as voltage and current 
waves (instead of electromagnetic waves), using elemental circuit theory. 

From (5.96), according to the fundamental property of the gradient, in 
the transverse plane the line integral of the electric held is path independent 
and consequently voltage V and potential difference <f> 2 — 5>i will be the same. 
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Figure 5.2: two conductor transmission linecambiar ejesj \ Comprobar si los ejes 
y el texto coincide. Los conductores se ponen en negro. entero no^ decir que 
cl son los conductores y que los sentidos de la corrinte in idreccion opuesta en 
cada conductor son indicadas. Poner origen coincidiendo con el conductor] jver 
dibujos de Salva. en todo ccaso la flecha de I de acuerdo con libro de siempre 
es al reves 



Then for TEM waves (using, without loss of generality, Cartesian coordinates) 
we have 

V = $(2) - $(1) = - / E t .dl = - ( E x dx + E y dy (5.99) 
Ji Ji 

where $(2) and <&(1) are the values of the scalar function <E> at the the con- 
ductors 1 and 2 and where I is any line that joins the equipotential transverse 
sections of these conductors (fig. 5. 2). Deriving with respect to z and taking into 
account Faraday's law, (1.1c), particularized for the source-free region, outside 
the conductors, we get 



dV 

dz 



-^-dx+-^-dy = --J-B v dx + B x dy (5.100) 



Note that J l —B y dx + B x dy is the magnetic flux through the area swept, 
along a unit of length in the direction z, by the line I joining the conducting 
surfaces. This flux can be expressed by using the magnetostatic definition of 
coefficient L of self-inductance per unit of length, as the product LI. Therefore 
we have 

dV dl 

IS = ~ L ai < 5 - 101 ' 

On the other hand, from Ampere's law, (l.ld), for the source-free dielectric 
region, we have 

' H t .dl= (f H x dx + H v dy (5.102) 
i Jr 
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VW — orrY^r 




Figure 5.3: 



where T is a closed path around one of the wires (see Fig. 5.2). Deriving with 
respect to z, we have 



w z = j^- dx+ ^ dy = -wti Dxdy ~ Dydx (5 ' 103) 



where, in a similar way as above, —D y dy + D x dx represents the flow of vector 
D per unit length in the direction z. Using the magnetostatic definition of 
capacitance C per unit length, this flux can be expressed as the product CV . 
Thus we have 

dl dV 

Hi = - c m < 5 - 104 > 

Note that, from (5.99) and (5.102), V and I must have the same z dependence 
as E and H, respectively. Thus, V and I are also traveling waves. 
In summary, according (5.101) and (5.104) we have 



d® _ t dl 
dz ~ ^ dt 



dl _ ndV 
dz — ° dt 



(5.105a) 
(5.105b) 



which are the coupled differential equations that voltage and current satisfy at 
any z cross section of an ideal line composed of perfect conductors separated 
by a lossless dielectric. Equations (5.105) are called ideal "transmission line 
equations". The use of these equations to study the propagation of TEM 
waves in transmission lines is considered in Chapter ??. 

5.2.4 Boundary conditions for TE and TM modes on per- 
fectly conducting walls 

For a guiding transmission system with perfectly conducting walls, the general 
boundary conditions on the walls require that the tangential component Ex of 
E and the normal component H n of H be null, i.e. 

E T = nxE=0 (5.106a) 
H n = n ■ H = (5.106b) 
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where h is the unit vector normal to the conducting walls. However, for TM 
and TE modes, as shown below, these conditions can be simplified and reduced 
to equivalent ones which are expressed only in terms of the z component of the 
fields. For example for TM modes, the requirement that 

E z =0 (5.107) 

on the perfectly conducting guide walls suffices to ensure that Eqs. (5.106) are 
fulfilled. From (5.86a) and the gradient properties, we can see that E t is normal 
to the lines where E z — cte and, therefore, to the boundary of the conductor, 
since it represents a line with E z = 0. Given that E t and H t are perpendicular 
to each other, the magnetic field is tangential to the conductor and thus E z = 
is equivalent to Eqs. (5.106). 

For TE modes, the necessary and sufficient condition to ensure that Eqs. 
(5.106) are fulfilled is that the normal derivative of H z be null on the perfect 
conducting parts of the guiding structure. That is 

BT-f 

—^ = \7H z -n=(\7 t + V z )H z -h = (5.108) 
on 

where we have divided V into its transverse and axial components. Taking into 
account (5.91a), we see that 

V t H z -h = H t -h = (5.109a) 

which means that H t is tangencial to the conductor and therefore, due to the 
perpendicularity of the fields, we have 

hxE t = (5.110) 

In summary, the necessary and sufficient boundary conditions on the perfect 
conducting walls of the propagation system are 

Boundary conditions on the perfect conducting walls 



For TM modes ,_ . , 

E z =0 < 5 - ma ) 

For TE modes ,„ _ _ . 

dH z _ n (5.111b) 

5.3 Cutoff frequency 

From (5.76) and (5.79) we see that, for propagation to exist, (3 g must be real, 
and consequently, 

k 2 > h 2 (5.112) 
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For this reason, /3 C , denned as 

I3 c = h= 2 -^ (5.113) 

A c 

is called the cutoff wavenumber, and A c is the cutoff wavelength. Thus, from 
Eq. (5.79), we have 

Pl=k 2 -fc (5.114) 



a 

and, consequently 



T2 = 72 - 72 (5-115) 
A g A A c 



where A is the wavelength of a plane wave in the unbounded lossless dielectric 
medium filling the ( / /better guiding structure// /) waveguide, and A s is that of 
the wave in the guide. Thus we have 



k= 2 i; /3 C = ^; = £ (5.116) 

A A c An 



The cutoff frequency f c is defined 2 as 

/c - 2n ~ 2n^JIs ~ 2tt [bAU) 

where v p = u/k is the phase velocity in the unbounded medium filling the ( 
ijbetter guiding structure; ij) waveguide. Thus, from (5.79), the wavenumber /3 
can be expressed in terms of f c , as 



/? s = fcyi-(y) (5.H8) 

and the corresponding wavelength A 9 in the ( abetter guiding structure} //) 
guide is 

2tt A 

Jo 



\ = -= , (5.119) 



s 



JJt)' 



which is greater than A. According to (5.118) the wavenumber is imaginary for 
modes with frequencies below the cutoff frequency / c , i.e. / < f c (or A > A c ). 

These modes, called evanescent modes, are attenuated and cannot propagate along the guide. Thus, 
( jjbetter guiding structurej j j) waveguides behave as high-pass filters for the TE and TM modes since 
they cannot transmit any of these modes for which the wavelengths, in the unbounded medium filling 
the ( jjbetter guiding structurej jj) waveguide, exceed the value of the cutoff wavelength. 



2 For a guiding transmission system with more than one dielectric the cutoff frequency can be 
defined in a different manner than (5.117). See for example Section ??. 
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In terms of the cutoff frequency, the expressions of the wave impedances for 
the TM and TE modes (5.88) and (5.93) for Ztm and Zte become, respectively 



From (5.120a) and (5.120b), we can see that Ztm < V and Zte > V and they 
become imaginary below the cutoff frequency. Thus, for / < f c , the waveguide 
behaves, in this respect, as a reactive impedance. 
From (5.79), we obtain the dispersion relation 



which is analogous to that obtained in (??) for the transverse electromagnetic 
waves in a nonmagnetized plasma. The plot of the phase constant as a func- 
tion of the frequency to (dispersion diagram) is shown in Figure 5.2.3. The 
transversal broken line corresponds to oj c = 0, i.e. to an unbounded lossless, 
nondispersive medium in which the wave propagates at the phase velocity v p 
regardless of its frequency. The solid-line curve represents Eq. (5.121) and 
shows that the waveguide is very dispersive close to the cutoff frequency uj c . For 
frequencies lj >> uj c such that their wavelengths are much smaller than the 
transversal ( abetter guiding structure] //) waveguide dimensions, the walls do 
not affect the propagation and the velocity tends to v p . 




JV: texto Dispersion diagram. Es la figura de plasmas, ver tb pp 444 del Jonk 
The group velocity, v gg , within the guide is given by 




(5.120b) 



(5.120a) 




(5.121) 




(5.122) 
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which is smaller than the phase velocity v p in the unbounded medium. The 
phase velocity within the waveguide( / jbetter guiding structure// /), v pgi is 
given by 




(5.123) 



which is always higher than that in the unbounded medium and is frequency 
dependent. Hence single conductor ( abetter guiding structure///) waveguides 
are dispersive transmission systems. Note that 

V P9 ■ v gg = v l (5-124) 

For TEM modes, from (5.94), we have j3 g = k which is real and independent 
of the frequency. Thus, all frequencies propagate along a lossless transmission 
line at the same phase velocity v p as that of the unbounded homogeneous di- 
electric filling the waveguide and there is no cutoff frequency. 

5.4 Attenuation in guiding structures 

For a propagating mode an attenuation constant a, Owing to energy dissipation within the 

waveguide, can arise from losses in the non-perfect conducting walls (a c ) and 
in the non-perfect dielectric filling the waveguide (ad). Thus, the attenuation 
constant a consists of two parts a = ad + a c . Dielectric losses are generally 
negligible when ( abetter guiding structure! //) waveguides are filled with air, 
which has a lower dielectric loss than do conventional dielectrics. 
hay que decir TE y TM y TEM.. 

First, we analyze the losses for TE and TM due to a non-perfect dielectric 
and afterward the ones due to non-perfect walls. In any case, as generally occurs 
in practice, these losses are assumed to be very small. 

5.4.1 TE and TM modes. 

Dielectric Losses 

When the dielectric filling the waveguide is lossy the attenuation can be easily 
taken into account if in the expressions obtained for ideal dielectrics the real 
propagation constants k and f3 g are replaced by — j-f and — jj g , respectively, 
where 7 = a + j/3 and j g = a d + j(3 g are the complex propagation constants 
in the unbounded dielectric filling the waveguide and in the waveguide, respec- 
tively. Then, from equations (3.17), (5.79) and (5.113), we have 



fc 2 (l-jtan^) = - 7 2 = -7'+/3c 



(5.125) 
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Using the above expressions for 7 and j g and neglecting the term a 2 d , because the 
attenuation constant a d is very small, we find 

Pi = k 2 -p 2 g (5.126a) 

k 2 P 2 + 01 

a d = ^"tanJd = - JL — — -t&nSd (5.126b) 

z Pg l Pg 

Thus, the attenuation factor is proportional to the loss tangent, tan 8 d , of the dielectric fill- 
ing the waveguide. On the other hand, (5.126a) coincides with Eq. (5.114) for 
waveguides with ideal dielectric, and consequently the phase constant (and thus 
the wavelength) remains practically the same as those for a lossless waveguide. 
The dependence of the attenuation factor a d on the frequency (assuming a range 
of frequencies in which the permittivity of the dielectric remain unchanged) can 
be deduced by substituting the expressions of t&nSd and (3 g , given by (1.91) 
and (5.118), respectively, in (5.126b). Thus we get 



ad 2(l-(fjff) 1 ' 2 



(5.127) 



where r\ is the intrinsic impedance of the dielectric given in (3.33a) and a e is its 
effective or equivalent conductivity (f .78). From (5.127) we can see that ay becomes very 
high at frequencies close to the cutoff value, then decreases to a minimum value, 
and afterwards increases with the frequency, becoming almost proportional to 
it. 



Wall losses 

When the conductivity is finite the tangential magnetic field induces currents 
which are not restricted to the surface and, according to Ohm's law, are asso- 
ciated with a tangential electric field (i.e. J = aE = nx H) in the walls. The 
vector product of the fields E and H at the surface of the walls represents a 
flux of power directed towards the inner of the wall. This power coincides with 
the dissipation in the conductor caused by the Joule effect and is subtracted 
from the mode that propagates along the waveguide. As a consequence, the 
amplitude of the electric and magnetic fields of the mode are attenuated accord- 
ing to e~ a " z , where a c is the attenuation constant due to wall losses. We can 
determine the value of a c for a given propagating mode by taking into account 
that the time-average power P av transmitted through the cross-section S of the 
guiding transmission system is 

Pav = [ Vav-ds = l f Re(E t x H* t ) ■ ds (5.128) 
Js 1 Js 

Since, due to the losses, the amplitude of the field wave varies according to 
e~ a " z , then, P av will vary according to e~ 2a " z . Moreover, the law of energy 
conservation requires that the rate of the decrease of P av with distance along 
the transmission system equals the time-average power loss on the surface of the 
walls per unit length, P' d , in the direction of propagation. Therefore, we have 
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dz 



and thus 



a c = 



P' 

^d 



(5.129) 



(5.130) 



If H is the ^magnetic?? field existing near the walls, the time-average power 
dissipated per unit of length in the walls is given, according to (3.45), by 



H H dl 



(5.131) 



where R s is the surface resistance given by (3.46) and T is the cross-sectional 
contour of the non-perfect conducting walls. Thus the coefficient of attenuation 
of the n-th TE or TM mode is found to be 



Or: 



R s f r HH*dl_ R s J r HH*dl 
4 f s Pavds ~ 2 f s Rc(E t xH* t )-ds 



(5.132) 



This equation will be applied in next chapter to the calculation losses for TE and TM modes in 

waveguides. In strict terms, the modes we have found assuming perfect conducting 
walls are no longer valid since non-perfect conducting walls represent a change 
in the boundary conditions because in this case the tangential component of 
the electric field is not null. However, if the losses are small, we can make an 

approximate analysis (known in Mathematical Physics as "first order perturbation method") 

by assuming that the field configurations or modes in the waveguide coincide 
with those found for ideal-wall ( abetter guiding structure/ //) waveguides. 



5.4.2 TEM modes 

The coupled differential equations (5.105) for ideal transmission lines can be easily extended to lines 
with a non perfect of dielectric (constitutive parameters ^i, (7) separating the perfect conductors. 

In this case, at any z cross-section of the line, an additional current increment 
AI leaves .. 

assuming that the dielectric has a conductivity o such that 

AI = gV (5.133) 

in which g deno tes the conductance 

pp 487 del Jonk: the series and shunt low-parameters r and g... 



<9$ 

~dz~ 

dz 



L di 
dt y 



(5.134a) 
(5.134b) 
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Chapter 6 



Some types of waveguides 
and transmission lines 

6.1 Introduction 

In the previous chapter we examined some general properties of the propagation modes that may 
exist in an ideal guiding transmission system which has no sources and is constituted by perfect 
conductors and one ideal homogeneous dielectric. Specific expressions for such modes can be deter- 
mined only when the particular geometry of the guide is given. In this chapter we will first analyze 
in some detail the homogeneously filled rectangular and circular metallic waveguides. After this, 
as a simple example of non homogeneous guiding structure in which the electromagnetic field prop- 
agates in more than one dielectric, we will study the dielectric slab waveguide. Then, we will give 
some basic ideas on propagation in strip and microstrip lines. Finally, we will consider cavity res- 
onators which are basically constituted by a dielectric region totally enclosed by conducting walls. 
This region, when excited by an electromagnetic field, presents resonance with a very high-quality 
factor Q. In particular, we will study the common simple cases of rectangular and circular cavity 
resonators. 

6.2 Rectangular waveguide 

Figure 6.1 shows a rectangular waveguide of sides O, and 6, with d !> b. and homogeneously 
filled with a perfect dielectric. Following the theory developed in the previous chapter, in order to 
calculate the TE and TM modes that can propagate in this waveguide, we start by solving the wave 
equation for the longitudinal components Z of the field with the corresponding boundary conditions 
determined by the geometry of the system. The transverse components are then calculated from 
these longitudinal ones. 

With axis chosen as shown in the figure, the expressions for the fields in (5.76), take the form 

E = E {x,y)e- 3 ^ z (6.135a) 
H = H {) (x,y)e- j ^ z (6.135b) 
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Next, we are going to find the expression for these fields, first for TM modes and afterwards for the 
TE modes. 



6.2.1 TM modes in rectangular waveguides 

For the TM modes, the differential equation (5.80) for 

E z = E 0z (x,y)e~ 3l3 n z (6.136) 

Can be Solved by USing the standard method of separation of variables in rectangular 

Cartesian coordinates. For this, we assume, for Eq z , solutions in the form of the 
product 

E 0z (x,y) = X(x)Y(y) (6.137) 

in which X(x) and Y(y) are, respectively, functions only of x and y. 
By substituting (6.137) in (5.80) and dividing by Eq z , we get 

1 d 2 X 1 d 2 Y , 2 , 

h 2 = (6.138) 



X dx 2 Y dy 2 

As each summand depends on a different variable, it should be verified that 

' " ? " A -hi (6.139a) 

-h 2 y (6.139b) 



X dx 2 
1 d 2 Y 



Y dy 2 

h 2 = Pl = h£ + h 2 y (6.139c) 

where we have substituted, according to (5.113), h by the cutoff wavenumber 
f3 c and where h x and h y are the separation constants to be determined from 
the boundary condition (5.111a) at the guide walls. This boundary condition 



6.2. RECTANGULAR WAVEGUIDE 
for the geometry of Figure 6.1 implies 



119 



E 0z = at 



(6.140) 



y = 



The solution of the Eqs. (6.139a) and (6.139b) are, respectively, 

X = C\ s'm.h x x + C<i cosh x x (6.141a) 
Y = Cs sin h y y + C4 cos h y y (6.141b) 

where the coefficients are arbitrary constants to be determined from bound- 
ary conditions. Therefore the general solution (6.137) for E 0z takes the form 

Eoz = (Ci sinh x x + C2 cos h x x) (C3 sin h y y + C4 cos h y y) (6.142) 
From the boundary conditions (6.140), we find that C2 = C4 = and 



h x = 



7rm 
a 

K = T 



(6.143a) 
(6.143b) 



and thus, from (6.139c), 



(6.144) 



where 771 and 71 are integers. The different soiutions achieved by giving values to 771 and 71 are 
termed TM mn modes and each set of values of 771 and 71 indicates a specific mode. Thus, from 
(6.136), (6.142) and (6.143), for TM mn modes, we have 



TlTTl 



7T71 



E z = A mn e -"Vsm xsm—y 

a b 



(6.145) 



where the product of the constants C\ and C3 has been replaced by a new 
constant A mn . 

Once we know the longitudinal component E z , we can calculate the trans- 
verse components E t by means of (5.86a) and then, by using (5.86b), which 
implies that 

Z TM = ^ = -^ (6 .i46) 



H, 



H, 



we can obtain H t . As a result, we get the following general expressions for the 
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TM TO „ modes in rectangular waveguides 
(E z ) TMmn = A mn e-^ sin ^xsin ?fy 

(E t ) TMmn = (-jA mn ^^e-^s z cos^xsinify) x- 
(jA mn ^e-^° z sin ^xcos ^yj y 

(H t ) TMra „ = {jA^n-'^e-^^ sin^xcos^fy) x- 
(jA^-^^e-^^cos^xsin^y) y 



6.2.2 TE modes in rectangular waveguides 



To analyze the TE modes, we can follow a procedure similar to that used for 
the TM modes but now solving for H z and imposing the boundary condition 
(5.111b), dH z /dn = 0, on the guide walls. This, for the geometry of Figure 
6.1, implies that 



dH z 

dx 
dH z 

dy 



= at 



at 



X 


= 


X 


= a 


y 


= 


y 


= b 



(6.148) 



Then, using (5.91a) and (5.91b), and after steps analogous to those followed for 
TM modes, we get 



TE m „ modes in rectangular waveguides 

(-H" z )TE m 



B mn t jl3gZ cos ™xcos ^y 



(H t ) TEm „ = (jB mn ^^e-^ z sin^xcos^fyj x+ 



COS 



— a; sin -5-j/ y 



(jB. 

(E t )TE 

(jB mnV ^?fe-^* z sin ^xcos ?f 



) 



jB mn rifi c J fe 



k.me-iPgZ COS Sin 21Ly 



) 



y 



3.149) 



Note that, for both, TM mfl and TE mn modes, the subindexes 771 and Tl, indicate the num- 
ber of half-wave variations of the field in the x and y directions, respectively. 
For a TM m „ mode with morn equal to zero, from (6.145), we have (E z )te 00 = 
and consequently, from Eqs (6.147), (Et)TE 00 = and (H t )xE 00 = 0- Hence, 
there is no TM mode in which m or n is equal to zero. This was to be expected 
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because a TM wave with E z = would degenerate to become a TEM wave 
which, as we saw in Subsection 5.2.3, cannot propagate within a waveguide. 

For TE mn modes, it is easy to see from (6.149) that either morn may be 
equal to zero but not both at the same time, since in this case the expression 
of (H z ) TEmn in (6.149) reduces to 

(H z )te„ = B 00 e-^* z (6.150) 

while (H t ) TEoo = and E = (E t ) TEoa = 0, such that only (H z ) TEgo exists. 
This field does not fulfil Maxwell's equations, since a time-varying field H 
should generate an electric field E. Therefore the TEqo mode cannot exist. 



Cutoff frequencies in a rectangular waveguide 

From (5.117), (6.139c), and (6.143), we see that the cutoff frequency for either 
a TE m „ or a TM m „ mode is given by 



(f) = 

\jcjmn 2 



(«) + U) 



(6.151) 



where v p is the phase propagation velocity of the wave in the unbounded 
medium filling the waveguide. The wavelength and wavenumber in the waveguide 
are given, respectively, by 



(A c ), 



(f) 2 +(f) 2 



/wry 

\~b) 



(6.152a) 



(6.152b) 



From (6.151) we see that the cutoff frequency of the modes depends on 
the dimensions of the cross-section of the waveguide. Values of the cutoff 
wavelengths and frequencies for several modes are 



(A c ) TEi 

(Ac)te 

(A C ) TE 



(A c ) 



TEi 



— 2a; (/c)tEi = 
= 26; (/ c ) TEqi 

= a i (/c)tE 20 = 



= (A c ) 



2a 

_ Vp_ 

= 2b 

Vp 

a 

2ab 



TM i 



(a 2 +b 2 ) i 



(fc) 



TEi 



(fc) 



TMi 



(6.153a) 

(6.153b) 

(6.153c) 

_ v p (a? + l?)* 
^ 2a6 

(6.153d) 



Note that if a = b the cutoff frequencies of TEi and TEoi and the two modes 
are equal except for a rotation of n/2. 
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The dominant TEio mode 



In practice, we usually wish to have only the mode which has the lowest cutoff 
frequency (^called fundamental or dominant mode??) propagating through 
the guide. Thus, in the case of a rectangular waveguide, if a > b, such that 
(/c)tei ^ (/c)te i> ^ ne waveguide is usually designed so that only the TEio 
mode can be propagated. The cutoff frequency of the dominant TEio mode is 
selected by means of the dimension a. The ratio of the cutoff frequency of each 
mode to that of the TEio mode as a function of a/b is plotted in Figure 6.2. We 
see that the separation of the cutoff frequencies for different modes is larger for 
higher values of the ratio of the o and b dimensions Note that if a ~ 26, then 
the cutoff frequencies of the modes TE i and TE20 are nearly the same and in 
the frequency range v/2a < / < v/2b only the TEio mode can be propagated. 
Moreover, if a > 2b, then (/ C ) TE20 < (/ c ) TEol ■ As we will see in the next Section, 
losses due to non-perfectly conducting walls increase as b decreases. Thus, to 
have the greatest frequency range in which only the TEio mode can propagate 
and, at the same time, to have the smallest losses possible, we usually choose 
the dimensions of the guide such that a c± 2b. Under this condition, only TEio 
modes will propagate in the frequency range (/ c )tEi < / < 2(/ c )tEi - F° r the 
dominant TEio mode the general expressions (6.149) simplify to those given in 
(6.154) where the constant Bio has been replaced by Hq. 
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Rectangular TE 10 mode 



5 = - 

a 

, _ Vp 

J c ~ 2a 



/W p -(f) 2 

H z = H o e- j0 o z cos lx 
H x =jH Q ^e-^o z smlx 

H y = 

E z =0 

E x = 
E y = -jHo^e-^9 z S m^x 



(6.154) 



6.2.3 Attenuation in rectangular waveguides 

Losses due to a non-perfect dielectric filling the waveguide and to non-perfect 
conducting walls can be calculated using the expressions (5.127) and (5.132), 
respectively. For a given mode, to obtain the attenuation due to dielectric 
losses, we simply need to use, in the formula (5.127), the value of the cutoff 
frequency of the mode, given by (6.151), and the values of the constitutive 
parameters of the dielectric at the work frequency. However, to find the the 
attenuation constant a c due to wall losses for any TE or TM mode, though 
not complicated, is quite laborious. Here, to illustrate the procedure, we will 
consider the particular case of the dominant TEi mode 



Attenuation of the TEi mode For the TEio mode, the integrals of the 
formula (5.132) can be calculated from the general expressions for the field 
components (6.154). Thus, for the denominator, we have 



TEi 



jjP*v)T:*x -dS=~J (E y H* x ) TEio dxdy 



( aH \ 



2 



LU[iab(3 (6.155) 



V 2tt J 

Regarding the integral of the numerator in (5.132), because in the dominant 
mode TEio in a rectangular waveguide the magnetic field H has only H x and 
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H z components, this integral takes the form 



jf H ■ Hdl = 2 S^£ (\H X \ 2 + \H Z \ 2 ) dx + £ (\H X \ 2 + \H Z \ 2 ) dy 



Using the expressions of H x and H z given in (6.154) and by operating, we obtain 



(6.156) 



H H dl = 2Hl 



2H 2 



2 V/c 



+ b 



(6.157) 



where the last expression is obtained from (5.118). By substituting (6.155) and 
(6.157) in (5.132) and after operating, we finally obtain the following expression 
for the attenuation factor (o: c )tEio 



(«c) 



TEi 



(l + f(t) 2 



'1<L X 2 * 



1 + 



Np/m 



(6.158) 

Following a similar analysis, we can show that the general expresions for the 
attenuation constant et c due to wall losses for any TE mn mode is 



(«c) 



TE r , 



where 

While for TM m „ mode is 




-^m 2 +n 2 



(6.159) 
(6.160) 

(6.161) 



These expressions show the dependence of the attenuation on the frequency. Reedactar: I jjCom- 
puted values of Ol c for a few TE m71 and TM mn modes are given in Figure 6.3In practice, surfaces 
imperfections, the value of Ol c may be greater than the theoretical values. This effect can be reduced 
using well polished walls. 
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Figure 6.3: Atenuacion en guias rectangulares: a commom characteristic It 
tends to infinite when / is close to the cutoff frequency, decreases toward an 
optimum frequency (minimum value of a c ) an then increases almost linearly 
with / 



